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ABSTRACT 


This  report  is  concerned  with  the  study  of  instabilities, 
generally  referred  to  as  POGO,  resulting  from  coupling  between 
liquid-fueled  propulsion  systems  and  structural  dynamics.  The 
report  deals  specifically  with  analysis  of  active  *nd  passive 
methods  of  sta.ility  compensation,  methods  of  predicting  POGO 
lj.nu.it  ojcac  amplitudes,  and  statistical  methods  of  s*r-‘  .atir" 
POGO  stability  margins.  Zn  addition,  experimental  '■'tudiea  were 
conducted  to  evaluate  certain  passive  compensation  methods. 
Specifically,  sealed  bag  accumulator,  gas  injection,  and  cavi- 
tating  venturi  concepts  were  studied. 
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SECTION 


SUMMARY 


1.  PROGRAM  DESCRIPTION 

The  research  work  effort  concentrated  on  analyzing  and 
evaluating  hitherto  untried  methods  of  POGO  stability  compen¬ 
sation,  developing  a  unique  method  for  predicting  POGO  stabil¬ 
ity,  and  developing  a  modal  for  explaining  the  existence  of 
stable  limit  cycles  observed  on  all  linearly  unstable  config¬ 
urations.  In  addition,  analysis  was  performed  on  three  con¬ 
cepts  of  POGO  compensation  (sealed  bag  accumulator,  gas  injec¬ 
tion,  and  cavitating  venturi)  to  establish  the  relationship  of 
dynamic  characteristics  to  specific  design  parameters.  Two  of 
these  concepts  were  further  studied  experimentally  to  evaluate 
the  analytical  conclusions.  Results  of  the  analyses  and  exper¬ 
imentation  are  briefly  summarized  below. 

2.  CONTROL  VALVE  COMPENSATION 

Stability  compensation  was  achieved  using  a  valve  to  mod¬ 
ulate  suction  line  flow  rates  in  direct  response  to  suction 
pressure  perturbations. 

The  study  revealed  that  suitable  choice  of  control  gains 
required  rather  precise  knowledge  of  both  suction  line  and 
structural  frequencies.  The  method  would  be  least  effective 
on  configurations  where  compensation  would  be  required  most; 
namely,  under  conditions  where  suction  line  and  structural 
frequencies  were  equal.  This  limitation  could  be  avoided  if 
the  control  valve  were  placed  in  the  diecmrge  lines  where 
phase  between  flow  rate  and  pressure  is  not  affected  by  system 
frequencies. 

3.  TURBINE  SPEED  CONTROL 

Studies  revealed  that  pump  discharge  flow  rate  could  be 
modulated  by  controlling  fuel  flow  to  t.ie  turbine  gas  genera¬ 
tor.  The  gas  generator/turbine  system  acts  as  an  amplifier 
whose  output  varies  the  main  pump  flow  rate  in  response  to  sig¬ 
nals  introduced  in  the  turbine  bootstrap  system.  Result®;  of 
the  analysis  showed  that  the  method  of  control  was  effective 
when  suction  line  frequencies  and  structural  frequencies  were 
sufficiently  separated.  However,  if  suction  l.ne  frequencies 


were  equal  or  nearly  equal  to  structural  frequencies,  this 
net hod  of  control  produces  equal  but  offsetting  effects  in  suction 
and  discharge  systems.  Thus,  under  the  conditions  of  equal 
suction  line  and  structural  frequencies  when  the  system  stability 
is  aost  critical,  the  Method  of  ccapensatioa  is  least  effective, 

k.  AMPUTEES  LIMITING  STUDIES 

Analysis  showed  that  stable  liait  cycles  in  the  amplitude 
range  observed  during  Titan  II  flights  could  be  simulated  by 
a  Mathematical  aodel  in  which  the  only  non-linearity  present 
arose  from  thermodynamic  processes  in  the  pump  cavitation 
region.  The  non-linearity  was,  however,  too  weak  to  make  liait 
cycle  amplitudes  independent  of  other  system  parameter*.  In 
fact,  the  study  showed  that  liait  cycle  amplitude  was  extremely 
sensitive  to  uncertainties  in  pump  gain,  structural  frequency, 
and  auction  line  frequency.  Normal  but  uncontrollable  variation 
in  these  quantities  from  one  vehicle  to  another  could  product 
an  order  of  magnitude  uncertainty  in  liait  cycle  amplitude 
prediction. 

5.  STATISTICAL  STABILITY  STUDIES 

A  method  was  developed  for  relating  system  gain  margin  to 
probability  of  syatea  stability.  Statistical  distributions  of 
critical  system  parameters  treated  as  random  or  functionally 
related  variables  were  combined  by  Monte  Carlo  analysis  tech¬ 
nique  which  simulated  any  specified  number  of  flights.  Results 
from  analysis  of  a  representative  Titan  III  MOL  configuration 
indicate  that  the  customary  method  of  root  sub  squaring 
tolerance  response  variations  results  in  gain  margin  require¬ 
ments  that  can  be  too  large  by  a  factor  of  3  or  4. 

6.  FLUIDIC  AMPLIFIER  COrtSGL 

On ry  preliminary  analysis  was  require  .  o  demonstrate 
that  fluidic  amplifiers  were  not  practical  dav^'es  for  modu¬ 
lating  main  stream  flow  in  response  to  pressure  or  accelera¬ 
tion  signals.  The  principle  disadvantage  of  the  fluidic 
amplifier  results  from  the  pressure  loss  that  must  occur  by 
virtue  of  requiring  a  rather  large  part  of  the  total  flow  for 
use  as  a  control  jet.  Also  pressure  noise  resulting  from 
■suring  the  control  jet  with  main  flow  can  mask  the  signal 
present  in  the  main  flow.  Because  of  the  pressui e  losse-, 
the  method  would  impose  intolerab1  -  performance  losses.  TV.? 
noiae  problem  would  make  the  concrol  insen8.it ire  to  pressure 


7.  SEALED  BAG  ACCUMULATOR 

The  sealed  bag  accumulator  POGO  fix  was  evaluated  to  deter¬ 
mine  the  effects  of  accumulator  design  parameters  (geometry  and 
syutem  static  pressure)  on  its  dynamic  characteristics  (iner- 
tance,  compliance,  and  resistance).  One  series  of  tests  deter¬ 
mine  that  accumulator  compliance  was  independent  of  geometry 
for  the  designs  tested  and  primarily  a  function  of  pressure. 

As  system  static  pressure  was  increased  accumulator  compliance 
decreased.  Through  dynamic  testing  the  mathematical  model 
synthesis  for  the  accumulator  was  found  to  be  in  error.  This 
result  inhibited  data  analysis  and  only  qualitative  results 
were  attainable. 

8.  STJv.iON  SYSTEM  GAS  INJECTION 

The  gas  injection  program  consisted  of  both  analytical  and 
experimental  phases.  Analysis  resulted  in  a  math  model  for 
predicting  suction  system  response  with  gas  injection.  The  ex¬ 
perimental  effort  consisted  of  five  test  sequences  utilizing 
water,  aerozine  50  and  nitrogen  tetroxide  as  propellants.  The 
experimental  results  indicate  that  a  correlation  between  anal¬ 
ysis  and  experiment  is  in  general  difficult  to  obtain  with  a 
horizontal  test  section.  Vertical  testing  yields  the  best  re¬ 
sults.  Additional  conclusions  are  that  certain  gases  -  CO2  in 
water  and  N2  in  nitrogen  tetroxide  are  dissolved  quite  rapidly 
and  produce  no  significant  shift  in  suction  system  frequencies. 
Furthermore,  it  was  found  that  gas  injection  does  affect  pump 
performance.  An  evaluation  of  the  extent  to  which  altered  per¬ 
formance  is  undesirable  must  be  considered  in  terms  of  the  par¬ 
ticular  turbo pump  application. 

9.  CAVITATING  VENTURI 

The  cavitating  venturi  which  through  previous  analyses  was 
found  to  induce  unacceptable  pressure  losses  in  the  suction 
system  was  re-examined  with  a  modification  which  it  was  felt 
would  reduce  losses.  The  modification  provided  for  a  low  loss 
bypass  tube  for  most  of  the  flow.  Through  analysis  it  was 
found  that  even  with  a  bypass  the  venturi  losses  were  still 
sufficient  to  make  the  device  unacceptable  as  a  POGO  fix. 


SECTION  II 


INTRODUCTION 


U  OBJECTIVES 

ill  program  objectives  sought  to  extend  and  improve  the 
basic  understanding  of  the  system  coupled  POGO  instability 
problem.  In  particular,  analysis  was  directed  toward  the  fol¬ 
lowing  ends  s 

&,  Further  refinement  of  methods  for  evaluating  and 
determining  system  stability. 

b.  Evaluation  and  comparison  of  relative  effectiveness 
of  passive  and  active  POGO  compensation  methods. 

o.  Refinement  of  models  for  critical  system  elements 
such  as  the  pump  and  auction  line. 

The  experimental  efforts  were  designed  to  complement  cer¬ 
tain  facets  of  the  analytical  studies  and  provide  means  of 
substantiating  analysis  conclusions.  The  studies  and  analysis 
that  were  made  a.r a  briefly  described  as  follows t 

1)  Turbine  Speed  Control  ~  feasibility  studies 
were  conducted  of  modulating  pump  flow  by  using 
auction  pressure  or  acceleration  signals  to  con¬ 
trol  fuel  delivered  to  the  gas  generator. 

2)  Active  Valve  Control  -  studies  we^-e  made  to 
determine  whether  POGO  stability  could  be  im¬ 
proved  by  controlling  suction  line  and  discharge 
flow  rats  in  response  to  suction  pressure  or 
acceleration  signals. 

5)  Fluidic  Amplifier  -  the  fluidic  amplifier  was 
investigated  to  determine  whether  suction  line 
or  discharge  line  flow  could  be  modulated  with¬ 
out  the  use  of  mechanics  valves, 

4)  Amplitude  Limiting'  Studies  -  analysis  was  per¬ 
formed  to  identify  the  non-linearity  responsible 
for  stable  limit  cycling  on  basical" v  unstable 
configurations. 


5)  Statistical  Stability  Analysis  -  a  method  was 
developed  for  basing  the  probability  of  system 
stability  on  various  critical  system  parameters 
treated  as  random  or  functionally  related 
statistical  variables. 

6)  Sealed  Bag  Accumulator  -  the  purpose  of  examin¬ 
ing  this  device  was  to  l)  determine  the  rela¬ 
tionship  of  accumulator  compliance  to  system 
static  pressure,  2)  devise  an  analytical  model 
which  would  relate  accumulator  dynamic  charac¬ 
teristics  to  specific  design  parameters,  and 

5)  verify  the  results  of  the  analytical  effort 
with  experimental  data. 

7)  Gas  Injection  POGO  Compensation  -  this  method  of 
correcting  POGO  was  studied  in  order  to  l)  de¬ 
velop  an  analytical  model  describing  the  suction 
system  performance  with  gas  injection,  2)  verify 
the  model  through  experiment,  3)  establish  the 
differences  between  horizontal  and  vertical 
testing  techniques,  4)  evaluate  the  effects  of 
gas  injection  n  pump  performance,  and  5)  estab¬ 
lish  the  effectiveness  of  certain  gases  for  POGO 
compensation. 

8)  Cavitating  Venturi  -  the  cavitating  venturi 
concept  was  re-examined  to  determine  if  a  by¬ 
pass  tube  would  reduce  pressure  losses  to  the 
extent  that  this  method  of  POGO  compensation 
would  become  feasible. 


SECTION  III 


ACTIVE  CONTROL  VALVE  ANALYSIS 


1.  SYSTEM  DESCRIPTION  AND  ANALYSIS  APPROACH 

The  introduction  of  Section  I  contained  a  discussion  of 
the  potential  advantages  associated  with  an  active  POGO  Control 
Device*  The  active  control  valve  concept  was  one  of  the  com¬ 
pensation  Methods  examined  during  the  Phase  I  Program.  Choice 
of  the  active  valve  as  a  method  of  compensation  originated 
from  recognising  that  it  would  be  desirable  to  control  directly 
the  flow  of  propellant  to  the  combustion  chamber.  All  analysis 
was  conducted  on  typical  Titan  III  vehicle  configuration*  in 
order  to  use  readily  available  parameter  values  and  to  assure 
that  results  would  not  be  subject  to  too  asny  idealizing 
assumptions  required  of  a  more  general  analysis.  A  wide 
enough  range  of  actual  parameter  valves  for  stricture  and 
propulsion  system  were  considered  to  assure  that  conclusions 
drawn  ware  not  based  on  a  special  situation  peculiar  to  any 
particular  set  of  parameter  values. 

a.  Analysis  Methods  -  The  first  part  of  the  analysis 
was  conducted  on  the  analog  computer  with  s  complete  simulation 
of  the  Titan  III  B  structure  and  Propulsion  Systems.  The  pur¬ 
pose  of  this  effort  was  to  evaluate  the  feasibility  of  con¬ 
trolling  POGO  with  active  butterfly  control  valves  located  in 
the  fuel  and  oxidizer  suction  lines.  In  addition  to  the  analog 
study  a  digital  analysis  was  performed  which  consisted  of  both 
frequency  response  and  root-locus  studies  of  the  decoupled  fuel 
end  oxidizer  systems. 

For  both  phases  of  Investigation,  valve  actuation 
was  cot-pled  to  sensed  suction  pressure  through  a  mechanical 
linkage  characterized  by  a  second  order  system  representative 
of  dynamics,  a  gain,  and  a  tine  lag.  For  a  detailed  description 
of  the  control  valve,  the  reader  is  referred  to  Appendix  1, 
Figure  1  shows  the  essential  blocks  of  the  oxidizer  side  of 
the  system  along  with  the  coupling  networks.  The  fuel  side  is 
similar  except  for  values  of  characterizing  parameters.  The 
control  valve  loop  is  Indicated  by  the  heavy  lines  In  the 
Figure.  The  followir  ■*  paragraph*'  contain  a  discussion  of  the 
math  model  development  along  with  the  details  and  results  of 
the  two  phases  of  analysis. 


Oxidizer  System  Block  Di»grau 
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MATHEMATICAL  DEVELOPMENT 


The  aquations  describing  the  dynamic  operation  of  a 
butterfly  valve  ware  derived  and  initially  combined  with  the 
discharge  line  equation  for  discharge  flew  modulation, 

a*  Discharge  Systeei  -  From  the  development  presented 
in  Appendix  II,  the  pressure  upstreaa  of  %  valve  (P^Q)  is 

related  to  the  downstream  pressure  (Pjq),  steadf  state  flow 

rate  (W  ),  steady  state  flow  coefficient  (C  ),  oscillatory 
V  •  ▼ 

flow  rate  (Wv<j),  *ad  oscillatory  flow  coefficient  (C^)  by: 


where  Cq  is  a  constant  allowing  for  different  fluid  specific 

gravities.  Froa  Appendix  I  the  equation  relating  oxidizer 
discharge  pressure,  flow  rate,  and  iapedance  is: 


Pj  -  P  >  L  ^  t  R .  W, 
do  c  d  do  d  do 


(2) 


Where:  Pdo  is  the  oscillatory  discharge  pressure 

Pc  is  the  oscillatory  combustion  chaaber  pressure 

I.  is  the  discharge  line  inertance 
CL 

is  the  discharge  system  resistance 

• 

W.  is  the  oscillatory  discharge  flow  rate 
ao 

Upon  ending  the  valve  pressure  drop  (Equation  1)  to  the  iner- 
tive  and  resistive  pressure  drops  of  the  discharge  system, 
Equation  2  becomes: 


The  above  relationship  then  describes  the  behavior  of  the  oxi¬ 
dizer  discharge  system  containing  an  active  valve  element.  A 
similar  relationship  exists  for  the  fuel  discharge  system. 


b.  Suction  System  -  Equation  1  was  subsequently  com¬ 
bined  with  the  suction  system  equation  as  described  in 
Appendix  II,  Repeated  here  the  resulting  functional  relation¬ 
ship  is: 


I  V 
so  so 


2f  co  I  V 

so  SO  80  80 


80 


' ^80*80 


ITo%o) 


-  p,o  *  [2  t.'Voho  b  *  2“v 
j\  +  «C(1  *  V  •  %  *  VJ 


(4) 


Where:  I  is  the  suction  line  inertance 

so 

r  is  the  suction  line  damping  ratio 
so 

m  is  the  suction  line  natural  frequency 
80 

• 

W  is  the  steady  state  suction  flow  rate 

is  the  steady  state  7alve  flow  coefficient 

L  is  the  suction  line  length 
so  & 

g  is  the  absolute  acceleration  of  the  suction  line 
so 

is  the  height  of  propellant  in  oxidizer  tank 


g_  is  the 
in  the 

P  is  the 
so 

P  is  the 
o 

A  is  the 
so 

is  the 

is  the 

«  • 

and  Qg  are  respectively  the  first  and  second 
mode  gimbal  velocity. 

C  is  the  oscillatory  valve  flow  coefficient  which 
is  a  function  of  blade  angle. 

c.  Valve  and  Suction  System  Characteristics  -  As  dis¬ 
cussed  in  Appendix  1 1  the  steady  state  valve  flow  coefficient 
(C  )  and  the  oscillatory  flow  coefficient  (C  )  will  be  a 

V  vo 

function  of  the  steady  state  butterfly  position.  Figure  2 
shows  the  variation  of  both  steady  state  and  oscillatory  flow 
coefficients  aa  a  function  of  blade  angle  for  a  typical  butter¬ 
fly  valve  used  in  this  investigation.  The  steady  state  blade 
angle  was  c’ncaan  as  60  degrees.  When  Equation  (4)  is  combined 
with  equations  of  Appendix  I  and  valve  control  loop  added,  tae 
suction  system  equation  upon  taking  the  Laplace  transform 
becomes: 


absolute  acceleration  of  the  propellant 
oxidizer  tank 

oscillatory  suction  pressure 
oxidizer  density 

cross  sectional  area  of  the  suction  line 
pump  velocity  relative  to  the  gimbal  point 
tank  bottom  modal  amplitude  ratio 


so 


f  M  S  +  C-- 
SO  80  1C 


*s) 


r  i  ,  i*1  1  s2 . 

L1..-..  (zdo  *  k=o  +  vsa 

*  Pao  -  Vto0#  (s2  t  2“„s  t  .72)("rd^  ToS^)  ‘  l5’ 

^80  +  S’o'1  +  S  Jso  “ao  p0^80  +  ^10  o^eoj^G  *  2 
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Where: 


is  the  steady  state  flow  coefficient 


W  is  the  steady  stst*  flow  rate 
vo  J 

*  2  °o  2 

C_  is  the  one  History  flow  coefficient  -rSM  ■%- 

w  V  v 

V 

w  is  the  natural  frequency  of  the  control  system 


f  is  the  damping  ratio  for  the  control  system 

is  the  effective  control  system  mass 

AQo  is  the  control  loop  gain 

Tq  is  the  control  loop  time  constant 

The  fuel  suction  system  equation  corresponding  to  Equation  (5) 
ia: 


'sf 


Mf  +  1 


if  *  ««  *  Kcf  *  V)SJ 


“sf5  *  2  Wsf5  *  V' 


+  P  -  c 
sf  If 


i2  )  (\{l  <6) 

V  V  V 


'A(l.f  *%>*l  <2  «f  “st  Vw  •  r +  c 


If 


f^f^G 


2g, 


The  fuel  ar  ncidizer  equation  defines  the  complete  relation¬ 
ship  betw  yu»p  suction  pressure,  P  ,p  ,  valve  flow  rate, 

•  *  SO  sf 

W  ,W  .,  discharge  system  parameters,  Z,,  t'.  ,  R  ,  and  para- 
vo  vi  a  c  p 

meters  affecting  the  control  dynamics,  AQo,  M^, 

The  equations  wsre  incorporated  in  one  complete  FOGC  system 
mathematical  model  end  programmed  for  stability  analysis  on 
an  analog  computer. 


1? 


3.  ANALOG  COMPUTER  ANAli'CIS 

As  a  first  step  of  the  analog  investigation,  it  was 
necessary  to  evaluate  the  behavior  of  the  system  without  the 
active  control  valves.  This  provided  data  which  would  serve 
as  a  basis  for  evaluating  control  system  behavior. 

a, .  Basic  No-Fix  Cystem  -  Appendix  I  presents  the  basic 

uncompensated  system  equations  which  were  programmed.  Figure 
3  displays  the  network  coupling  the  various  subsystems  to¬ 
gether.  For  evaluation  of  both  the  basic  system  and  the 
active  control  valve,  various  Till  configurations  were  chosen 
corresponding  to  flight  times  of  80,  120,  and  148  seconds. 
Figures  4,  5i  ari  6  illustrate  for  the  above  times  the  char¬ 
acteristics  of  the  basic  uncompensated  system  in  terms  of 
Nyquist  plots.  The  system  loop  has  been  opened  at  thrust  T, 
and  the  Nyquist  plots  are  for  the  gain  and  phase  of  ths 
quantity  T  ^/T^  (Reference  Figure  3).  Although  the  three 

figures  show  the  uncompensated  system  to  be  stable  (encircle¬ 
ment  of  the  +1  zero  Phase  gain  indicates  an  unstable  condition) 
the  overall  system  gains  are  sufficient  to  demonstrate  any 
stabilizing  effects  produced  with  the  control  valves. 

b.  System  with  Control  Loop  -  The  system  equations 
used  for  the  control  valve  simulation  were  those  given  in 
Appendix  I  with  the  exception  of  the  suction  line  equations. 

The  auction  line  equations  were  '  eplaced  with  Equations  (5) 
and  (6)  above.  Before  proceeding  with  evaluation  of  the  com¬ 
pensation  devices,  it  was  necessary  to  determine  values  for 
the  control  loop  gains,  AQq  and  A^,  which  would  produce 

stabilization.  To  accomplish  this  task,  the  oxidizer  and  fuel 
systems  were  considered  separately  without  the  influence  of 
pump  dynamics.  Referring  to  Appendix  I,  the  fuel  and  oxidizer 
systems  are  decoupled  by  replacing  Pc  with  PCQ  in  the  oxidizer 

equations  and  P  with  P  ,  in  the  fuel  equations.  Upon  taking 

C  Cl 

the  Laplace  transform  of  the  system  of  equations,  the  follow¬ 
ing  transfer  functions  car.  be  written. 
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OXIDIZER 

HJEL 

«*  *  V/*m 

<«*  *  ’4)/Tte 

(7a,  b) 

•  »  •  • 

V«4 1  V 

•  t  •  • 

fa/'5*  *  V 

(8a,  b) 

vp.« 

VPaf 

(9a,  b) 

p=Ao 

fA 

(10a,  b) 

?out/Pco 

T  „/P  , 
out,  cf 

(11a,  b) 

Combining  expressions  (?)  through  (11)  will  yield  the  open  loop 
transfer  functions, 

^out/^in'ox  ^out^in^fuei  ( L2a,  b) 

Substitution  of  the  appropriate  .flight  time  parameters  produces 
the  following  forms: 


OHDX2S3R  SYSTEM 

SO  Seconds:  w  =  69  Radians/Second 
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r* 
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120  Seconds:  w  -  8l  Radians/Second 
T 
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out 
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An  examination  of  Equations  (13)  through  (i>)  Indicates  that 

certain  values  of  both  positive  and  negative  control  gains  will 

reduce  the  oxidizer  T  ,/T.  fun^t  j  on**.  Scalier  positive  con - 

ou ;  in  r 

trol  gains  for  however  will  accomplish  a  greatei  total 

system  gain  reduction  than  will  negative  control  game.  Baaed 
on  these  results  oxidizer  Centro*  loop  gains  (A^;  of  +50  and 

+100  were  selected  for  the  initial  analog  anulysi*. 

Considering  the  fuel  system,  Equations  06;  through 
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(18)  again  indicate  overall  system  gain  reductions  can  be 
accomplished  with  either  positive  or  negative  control  gaine. 

In  this  case,  however,  negative  values  are  more  desirable 
than  positive.  Hence,  fuel  control  gaine  of  -50  and  -100 
were  selectsd  for  initial  examination. 

At  this  point,  the  question  should  be  raised  as 
to  what  system  characteristics  necessitate  gains  of  opposite 
sign  in  the  fuel  and  oxidizer  control  loops.  The  answer  is 
found  through  an  examination  of  the  sv  tion  system  pressure- 
flow  rate  relationships.  At  a  frequency  corresponding  to 
the  resonant  frequency  of  the  first  structural  mode,  the 
oxidizer  pressure  and  flow  rate  are  in  phase.  Thus,  with 
a  positive  control  gain,  an  increasing  pressure  drives  the 
valve  closed,  which  in  turn  reduces  the  amplitude  of  the 
flow  rate.  The  fuel  suction  pressure  and  flow  rate,  however, 
are  approximately  180  degrees  out  of  phase  at  the  structural 
first  mode  frequency.  Thus,  if  the  flow  rate  amplitude  is 
to  be  reduced,  a  negative  gain  is  required. 

The  dependence  of  control  gain  on  ratio  of  suction 
line  tc  structural  frequency  represents  an  undesirable  char¬ 
acteristic  of  the  control  system.  If  either  the  fuel  or 
tacidizer  suction  line  frequencies  crossed  over  the  structural 
frequency  during  flight,  a  change  in  sign  on  the  control  gain 
would  be  required.  Several  solutions  to  this  problem  exist 
and  wUl  be  discussed  in  later  paragraphs. 

4.  ANALOG  ANALYSIS  RESULTS 

The  analog  data  was  reduced  to  the*  form  of  Nyquist  plots 
which  are  shown  in  Figures  7  through  12. 

a.  Oxidizer  Fix  with  no  Time  Lags  -  Figures  7,  8,  and 
9  show  the  effect  of  an  oxidizer  valve  alone  on  system  response 
for  flight  times  of  80,  120,  and  '  ^8  seconds  respectively.  The 
control  loop  in  these  cases  is  laracterized  by  a  gain  (A ^  ) 

of  +50  and  zero  time  lag.  Compared  with  the  no-fix  Nyquist 
plots,  it  is  observed  that  the  oxidizer  valve  produces  signi¬ 
ficant  changes  in  gain  for  both  the  80  second  and  120  second 
cases.  Specifically  at  80  ssconda,  the  zero  phase  gain  is 
reduced  from  .27  to  .12;  and  at  120  seconds,  the  gain  is  reduced 
from  .55  to  .40.  At  148  seconds,  the  oxidizer  valve  produces 
little  effect  on  system  response.  Tnis  weak  influence  can  be 
explained  in  terms  of  the  reduced  coupling  between  the  oxidizer 
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system  and  structure.  Tne  strongest  coupling  at  148  seconds 
exists  between  the  fuel  suction  system  and  structure. 

b.  Oxidizer  Fix  with  Tise  Lags  -  Time  lags  added  t*  ths 
oxidizer  control  loop  alter  the  Nyquist  plots  to  those  shewn 
in  Figures  10  and  11  from  80  to  120  seconds  respectively. 

In  all  cases  the  time  lags  tend  to  reduce  the  gain  margin  im¬ 
provements  achieved  with  the  active  valve.  For  a  flight  time  of 
148  seconds,  time  lags  produced  little  effect  on  the  Nyqu:  st 
plots?  *g*i»*->  due  to  the  weak  coupling  of  oxidizer  and  structure 
systems. 

e.  Oxidizer  and  Fuel  Fix  -  Addition  of  the  fuel  control 
valve  to  the  system  produced  the  most  significant  changes  at 
120  and  148  seconds.  At  80  seconds,  the  coupling  between  the 
fuel  system  anu  structure  is  quite  weak;  hence,  as  in  the  case 
of  the  oxidizer  system  at  148  seconds,  the  fuel  oontrol  loop 
is  relatively  ineffective.  The  performance  of  the  fuel  fix 
can  be  assessed  from  the  Nyquist  plots  of  Figures  7,  8  and  12. 

d.  System  Instabilities  •*  During  the  analog  analysis,  a 
number  of  unexpected  results  were  encountered.  For  sxaaple, 
divergence  of  the  open  loop  system  response  occurred  when  an 
oxidizer  control  loop  gain  of  +100  was  used.  The  second 
problem  arose  while  attempting  to  evaluate  the  effect  of  time 
lags  in  the  fuel  control  loop.  In  this  case,  a  diverging 
oscillation  appeared  in  the  system.  Both  of  these  instabili¬ 
ties  can  generally  be  associated  with  an  unstable  internal 
system  loop.  The  only  csndidates  in  this  case  were  the  oontrol 
loops  themselves.  In  order  to  determine  the  particular  elements 
in  the  control  system  which  led  to  the  observed  instabilities 
a  root  locus  analysis  was  initiated  utilizing  an  existing 
digital  computer  program.  The  details  and  results  of  this 
analysis  are  presented  in  the  following  section: 

5.  DIGITAL  COMPUTER  ANALYSIS 

The  analysis  of  the  observed  control  loop  instabilities 
was  performed  using  the  Martin  Marietta  Corporation’s  Root 
Locus  Program  QD  030,  The  decoupled  fuel  and  oxidizer  P 

transfer  functions  provide  all  the  functional  relationships 
necessary  for  the  investigation. 

a.  Digital  Equations  -  Referring  to  Figure  13  below 
the  decoupled  P ^/g^  transfer  function  describing  the  suction 
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system  containing  an  active  control  vulva  is  expressed  as: 


V«b*  T-1 


where : 


a  Ph 
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In  the  above  expressions,  the  paras et era  have  been  written  in 
,a  general  fora  such  that  the  P /g^  function  applies  to  both 

the  fuel  and  oxidiner  systems. 
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Figure  13  Suction  Block  Diagram 
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Substitution  of  the  above  expressions  in  Equation  19  yields  a 
p  y'gg  transfer  function  containing  »  sixth  degree  MSM  poly¬ 
nomial  in  the  denominate.;.  The  stability  of  the  control  valve- 
suction  system  can  be  established  from  an  examination  of  the 
roots  associated  with  the  denominator  polynomial.  If  the  real 
parts  of  the  roots  all  lie  in  the  negative  half  of  the 
plane,  the  system  will  be  stable.  Boots  lying  in  the  positive 
half  of  the  MS"  plane  will  produce  either  a  pure  syatea 
divergence  or  an  oscillatory  instability  depending  on  whethfc* 
or  not  the  roots  contain  imaginary  components.  A  root  at  zero 
produces  a  metastable  system. 

b.  Root  Locus  Results  -  Figures  14,  15,  l6  and  17 

present  the  re«>ui<.-  digital  analyse  ir.  terms  of  root 

loci  with  variations  in  control  parameters.  Figure  14  demon¬ 
strates  the  behavior  of  the  oxidizer  side  a#  the  control  gain, 
A^,  varies  from  0  to  +100  0=0).  It  i  i  cbsvrved  that  the 

system  is  stable  up  to  a  gain  of  6o.5«  wetantabJe  at  60.5,  and 
unstable  for  A^  between  60.5  and  100  Alr.o  ft  will  be  noted 

that  the  roots  for  control  gains  of  6'  * !?  t.r  i  greater  have  no 
imaginary  component.  Hence,  the  inetabi.* will  produce  a 
pure  divergence  as  observed  during  the  analog  analysis. 

The  characteristics  of  the  oxidizer  control  loop 

are  shown  in  Figure  15  for  variation;;  in  control  loop  time 

lag  t  .  The  curves  indicate  the  nyste*  is  stable  for  varia- 
o 

tions  in  r  from  0  to  .050  (A^  «  +50){  however,  an  unstable 
condition  is  approached  as  r  is  increased  beyond  .050. 

Figure  16  which  demonstrates  the  behavior  of  the 
fuel  system  to  variations  in  AQj,  is  very  similar  to  Figure  14, 

For  the  fuel  valve  an  A ^  equal  to  84.5  defines  the  boundary 

between  a  stable  and  an  unstable  system.  Finally,  the  stability 
of  the  fuel  side  for  variations  in  tine  lag  is  observed  in 
Figure  17.  The  system  is  stable  for  1  o  0}  however,  as  r  in¬ 
creases  up  to  about  .059*  the  system  is  unstable.  In  this 
case,  due  to  the  lumginary  components  of  the  rocts,  the  system 
will  experience  an  oscillatory  instability.  Again  this  result 
is  confirmed  by  the  analog  analysis.  The  data  also  shows  that 
stability  is  regained  for  time  lags  between  .055  and  .050. 

c.  Valve  Effects  on  Suction  System  -  Jn  addition  to 
defining  regions  of  stability  for  the  system,  the  root  locus 
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Figure  15  ?J Poles  in  3"  Plan*  -  A,. ^  *  *>0 
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Analysis  demonstrates  the  manner  in  which  the  control  tsIts 
changes  the  suction  system  response.  In  Figure  14,  the  aagni- 
tude  of  the  iaaginary  root  coapoments  in  the  negative  half 
plane  is  a  direct  indication  of  suction  system  frequency. 

Thus,  it  is  observed  that  increasing  the  oxidizer  gain  froa 
aero  to  60.5  degrees  the  suction  systea  frequency  froa  9*3  to 
0.0  ops.  A  further  increase  in  introduces  a  negative  fre¬ 
quency  and  the  pure  divergence  which  is  indicated.  Also 
apparent  is  the  fact  that  suction  systea  frequency  is  quite 
sensitive  to  A^  for  A^  in  the  vicinity  of  oO.  On  the  fuel 

systea,  the  frequency  decreases  froa  18.65  cpa  for  b  0  to 

0«0  eps  for  •  84.3.  Again  the  sensitivity  of  frequency 

to  A^  increases  as  the  frequency  approaches  zero.  The  final 

observation  is  that  increasing  the  control  loop  gain  results 
in  an  increase  in  the  suction  systea  daspiag  as  indicated  by 
the  increasing  angle  between  the  j  axis  and  the  root  vector. 

It  any  be  concluded,  therefore,  that  gain  in  the 
control  valve  loop  produces  a  change  in  the  suction  systea 
frequency.  In  teras  of  the  total  systea,  increased  stability 
is  achieved  by  separation  of  the  suction  systea  and  the 
structural  frequencies. 

The  effiets  of  tine  lag  (r)  are  indicated  in 
Figures  15  and  17*  Tiae  lag  is  a een  to  affect  P JZq  gain  in 

both  the  fuel  and  oxidizer  systems.  In  the  oxidizer  side, 
increasing  r  initially  increases  the  systea  stability  up  to 
r  &  .015,  however,  additional  increases  tend  to  reduce  the 
stability  margin.  The  opposite  effects  are  observed  on  the 
fuel  systea.  As  t  is  increased  from  zero,  the  fuel  side 
becomes  more  unstable.  At  r  *  .007  *  maximum  in  instability 
is  indicated.  For  r  greater  than  .007,  the  roots  move  toward 
the  region  of  stability  and  the  auction  ayatem  control  loop 
ultimately  becomes  stable  for  r  greater  than  .035. 

6.  CONCLUSIONS 

The  purpose  of  this  investigation  was  to  determine  if 
HX30  Longitudinal  Oscillations  could  be  controlled  with  active 
valve  elements  in  the  vehicle  propulsion  systea. 
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a.  Control  Feasibility  -  Based  on  the  analyses  conducted, 
it  appears  that  this  net hod  of  control  is  in  fact  feasible.  The 
control  valve  affects  longitudinal  stability  by  changing  the 
suction  systes  natural  frequency.  The  frequency  shift  is 
acconplished  with  changes  in  the  control  loop  gain;  however, 
there  are  restrictions  inpcsed  on  the  gain  levels.  If  the 
control  gains  are  too  large,  a  divergence  in  the  systes 
characterised  by  the  control  valve  going  to  either  a  full  open 
or  full  closed  position  will  result.  The  sane  situation  will 
develop  in  the  fuel  systes  if  control  gains  are  too  large.  In 
addition,  the  fuel  system  appears  to  be  susceptible  to  oscilla¬ 
tory  instabilities  for  snail-tine  lags  in  the  control  loop. 

b.  Control  Disadvantages  -  One  of  the  disadvantages 
encountered  in  this  particular  control  system  is  ths  sensi¬ 
tivity  of  stability  to  control  valve  gains.  In  the  system 
studied,  it  was  necessary  under  sene  conditions  to  use  oantrol 
gains  of  opposite  sign  on  the  fuel  and  cod  diner  control  loops. 
This  requirement  w as  incurred  because  structural  first  soda 
frequency  are  higher  than  the  oxidizer  auction  system  frequency 
and  lower  than  the  fuel  system  frequency.  A  positive  gain  la 
the  oxidizer  control  loop  lowered  the  suction  frequency  thereby 
separating  it  from  the  structural  frequency  while  negative  gains 
in  the  fuel  system  achieved  frequency  separation  by  raising 

the  fuel  auction  frequency.  For  later  flight  times,  the 
structural  frequency  may  cross-over  the  fuel  auction  frequency 
thereby  requiring  a  change  in  aign  on  the  control  loop  gain. 
Furthermore,  it  appears  that  the  control  valve  will  have  little 
or  no  influence  in  the  system  when  suction  line  and  structural 
frequencies  are  equal.  Ther  *re  at  laaat  two  ways  of  avoiding 
the  above  difficulties.  If  «  parameter  which  maintained  the 
same  phase  relationship  with  flow  rate  were  used  to  actuate  the 
valve,  the  aign  change  could  be  avoided.  Such  a  relationship 
exists  between  the  suction  pressure  and  pump  discharge  flow. 

The  phase  inference  between  the  two  remains  at  approximately 
ldO  degree*  throughout  flight.  The  active  valve,  in  thia  caae, 
would  use  suction  pressure  for  actuation;  however,  it  would  be 
located  in  the  discharge  system  and  would  modulate  discharge 
flow  rate.  The  disadvantages  of  locating  the  valve  in  the 
discharge  system  as  presently  viewed  would  be: 

(1)  The  limited  physical  location  envelope  in  which 
the  valve  could  be  placed. 

(2)  The  possible  necessity  of  an  engine  modifica¬ 
tion  with  the  required  requalification  test  program. 


(3)  The  higher  pressure  losses  across  the  valve. 

The  other  possible  solution  is  to  use  control  loop 
gain  to  lower  both  fuel  and  axldixer  suction  ays ten  frequencies 
below  the  structural  first  node  frequency  for  all  flight  tines* 
One  probloa  which  would  be  encountered  here  is  the  eoctrene 
sensitivity  the  control  loop  frequency  displays  as  function  of 
gain  for  lew  frequencies.  Precise  knowledge  of  all  propulsion 
systen  parameters  would  be  required  in  order  to  establish  the 
proper  control  gain.  Unfortunately,  examination  of  the  latter 
FOOD  control  concepts  was  not  within  the  scope  of  the  present 
program.  If  the  active  control  valve  concept  is  to  be  pursued 
farther,  it  Is  rec amended  that  the  above  considerations  be 
evaluated  in  order  to  establish  an  optixun  control  valve 
configuration. 

f.  Suamary  -  In  conclusion,  the  active  valve  appears 
to  be  a  feasible  nethod  for  controlling  POQO.  There  are  a 
nunber  of  features,  however,  which  would  be  undesirable  for 
oertaia  vehicle  o  on  figurations.  It  is  felt  that  additional 
investigations  would  provide  practical  methods  of  avoiding 
these  disadvantages. 


SECTION  IV 


TURBINE  SPEED  CONTROL 


I.  GENERAL  THEORETICAL  CONSIDERATIONS 

a.  Introduction:  Control  of  POOO  Tin  Modulation  of  pro¬ 
pellant  flow  to  tho  turbine  gas  generator  represents  a  unique 
change  Iron  conventional  passive  methods  of  eliminating  coupled 
propulsion  system  flow  and  pressure  oscillations.  Control  is 
achieved  using  components  and  subsystems  already  present  in  a 
pump/ engine  assembly  and  adds  less  hardware  to  the  system  than 
most  of  the  other  passive  or  active  POGO  compensation  techniques. 
In  addition  to  providing  a  means  for  stabilizing  POGO,  turbine 
speed  control  offers  the  possibility  of  using  the  same  hardware 
components  to  vary  thrust  in  response  to  signals  generated  by 
trajectory  optimisation  guidance  commands. 

Control  of  the  system  is  achieved  by  using  the  turbinc/gas 
generator  as  an  amplifier  to  modulate  the  flow  to  the  main 
engine  combustion  chamber.  By  controlling  flow  to  the  turbine 
gas  generator,  turbine  speed  can  be  varied  which  in  turn 
directly  changes  pump  power  output  and  subsequently  p  tp  flow 
rate.  The  basic  cause  and  effect  relationships  follow  directly 
from  the  block  diagram  below. 


Figure  18  Basic  System  With  Turbine 
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As  shown,  the  control  elements,  Including  the  gas  generator 
supply  lines  (bootstrap  lines)  are  in  parallel  with  the  feed 
system  and  pump*  Control  is  obtained  by  driving  a  valve  in  the 
bootstrap  lines  with  a  signal,  g,  obtained  from  an  aeceleroaeter 
Control  cm  also  be  achieved  by  driving  the  gas  generator  valve 
with  a  suction  pressure  signal,  Ps.  Both  types,  acceleration 
and  suction  pressure  control,  will  be  discussed  subsequently* 

The  effectivity  of  either  method  depends  largely  upon  how  auch 
gain  attenuation  and  or  advantageous  phase  shift  can  be  obtained 
by  replacing  the  natural  feed  system  and  pump  dynamic  elements 
by  the  parallel  combination  of  those  elements  with  the  control 
dynamics  and  gas  generator/turbine  systems. 

b*  System  Model  Without  Compensation.  The  standard 
POGO  mathematical  model  does  not  ordinarily  contain  the  gas 
generator/turbine  dynamics*  The  impedance  of  the  bootstrap 
system  compares  with  the  main  discharge  system  impedance  is  too 
high  to  allow  POGO  induced  discharge  pressure  variations  to 
significantly  affect  the  bootstrap  systea  flow  rates*  That  is, 
referring  to  the  ricetch  below,  the  gain  around  the  loop  contain¬ 
ing  the  bootstrap  system  is  very  much  less  than  the  gain  around 
the  primary  loopt  hence  the  bootstrap  systea  is  usually  omitted 
from  POGO  analysis*  However,  if  the  gain  around  the  bootstrap 
system  loop  can  be  increased,  i*e*,  by  use  of  a  valve  that 
decreases  the  bootstrap  line  resistance,  then  that  feedbaok 
path  becomes  available  as  a  means  of  control* 


Figure  19  Basic  System  With  Control 


c.  Analysis  with  Gas  Gaaeratsr/Turbins  I*>op  Compensator. 

The  complete  set  of  equations  for  the  compensated  system  are 
developed  in  Appendix  ill  *  The  discussion  that  appears  below 
is  presented  with  the  aid  of  simplified  equations  and  block  dia¬ 
grams  that  retain  the  important  elements  of  the  system  and  control 
method* 


(l)  Control  via  acceleration  feedback  -  Since  It  is 
intended  to  stabilize  a  system  for  which  the  prime  consideration 
ia  limiting  oscillatory  acceleration,  acceleration  feedback  will 
first  be  attempted  as  the  priaar*7  control  gone*  The  block  diagram 
shows  the  pricipal  fuel  and  oxidizer  system  elements  including 
the  control  loop  through  the  gas  generator,  turbine,  and  pumps* 

As  shown,  control  is  achieved  by  direct  modulation  of  the  gaa 
generator  fuel  flow  by  the  acceleration  control  signal,  g.  Fuel 
flow  only  is  modulated  because  of  extremely  rich  mixture  ratio 
in  the  gas  generator.  The  transfer  function  from  g  through  the 
gas  generator,  turbine,  and  pump  is  given  by: 


aPdf*o  -  *c  aCr  °7  °4F.  a  „  5 

g  1  +  Q?  (\r  Qyjf  a8f  +Q4o  a3o  °8o) 


Where  the  elements  are  as  defined  on  the  block 
diagram  with  g  and  6 P^y,  o  the  incremental  change  in  fuel  or 

oxidizer  pump  discharge  pressure  per  unit  of  acceleration  input* 
When  the  transfer  function,  ^d  ,  is  reduced  to  essentials  in 

terms  of  Laplace  variable  S,  one  can  easily  observe  the  effect 
of  control.  Reduction  of  the  ^d  function  thus  produces  the 

following:  ® 


s 


K  G 


c 

1  + 


(a) 


Where  K  G  ia  the  total  control  path  gain  and 
®  S 

is  the  time  constant.  Such  a  simple  relationship  results  from 
the  fact  that  ia  much  larger  than  any  of  the  other  aystem 

time  o-  stants  so  that  they  can  be  neglected  for  purposes  of  il~ 
lustra  ng  the  control  method.  The  time  constant,  r  ,  ia  as- 

social  n  with  the  relatively  large  turbine  and  pump  rotary  inertias 
incorporated  in  the  system  element  G_.  Sffeeta  of  the  control  car. 
be  evaluated  by  observing  the  effect ^of  the  $^d  functions 

S 
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Speed  System  Complete  Block  Diagram 


in  the  main  eystea  functions  P  and  P  for  both  fu.l  and 

— SL  S— 

g  g 

oxidizer.  The  two  main  system  functions  can  be  simply  expressed 
if  the  main  fuel  and  oxidizer  systems  are  decoupled  at  F  to  per- 

nit  only  the  oxidizer  component  of  P_  to  feed  back  into  the  oxidizer 

system  and  only  the  fuel  component  of  P  to  feed  back  into  the  fuel 
system*  The  suction  system  response  with  control  present  is  then 
expressed  as: 


Where  0  is  the  total  control  function 
6 

P 

The  feed  system  response,  a  without  control 

g 

differs  from  the  response  with  control  only  by  the  presence  of  the 
control  term  in  the  numerator  of  the  transfer  function.  The  con¬ 
trol  term  also  affects  the  discharge  system  response  so  that  the 
total  propulsion  system  transfer  function  reflects  the  presence 
of  the  control  function.  Results  of  determining  the  effect  of 
control  on  the  combustion  chamber  pressure  response  expresses  the 
total  impact  of  the  control  function  as  follows! 


without  much  loss  of  generality  so  that, 
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K  /pH  - 
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_  Fro*  tht  transfer  function,  it  follows  that  the 
control  tern,  Qc»  cancels  under  the  condition  where  ■  O*  This 

condition  pertains  when  the  complex  frequency,  S,  is  equal  in 
Magnitude  to  the  suction  line  natural  frequency  since  Q  ■  s2  +  i. 

«2 

8 

Since  the  POGO  systen  instability  nanifests  itself  at  some  frequency 

close  to  the  natural  frequency  of  the  structure,  we  can  say  the  Q 

0 

■  0  when  the  natural  frequency  of  the  suction  systen  is  equal  tu 
the  natural  frequency  of  the  structure*  When  this  occurs,  the  Gq 

terns  in  the  P  transfer  functions  cancel.  This  is  equivalent  to 

JL 

f 

saying  that  the  control  would  not  be  effective  for  systems  where 
the  structural  and  station  systen  material  frequencies  were  equal 
or  nearly  equal. 


(2)  Compensation  for  Various  Configurations  Dependent 
on  Batio  of  Feed  Systen  Frequency  to  Structural  Frequency.  Although 
the  acceleration  control  method  is  not  effective  when  the  structure 
and  suction  line  natural  frequencies  are  equal,  system  stability 
can  be  considerably  improved  for  structural  nodes  with  frequencies 
different  froa  the  suction  line  modes. 


The  responses  will  be  discussed  in  terns  of  Nyquiet 
diagrams  or  frequency  response  plots  of  system  gain  and  phase  with 
frequency  as  a  parameter.  A  typical  plot  for  both  a  stable  and 
unstable  system  are  shown  in  Figure  21, 


The  maximum  response  vector  corresponds  to  a  fre¬ 
quency  that  is  very  nearly  equal  to  the  frequency  at  the  structural 
mode  one  is  attempting  to  stabilise. 


Figure  21  Typical  System  Response 


(3)  Structural  Frequency  Greater  than  Suction  Lina 
Frequency.  The  response  function  of  concern  is  - 


!*_  £2* .  !£2 .  i 
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For  purposes  of  discussion  the  total  open  loep 
systea  response  is  broken  into  fire  parts  designated  a  -  e,  where  - 
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The  response  is  then  given  by: 


“out 

*in 


*  ajb(c  +  d)  +  ej 
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Typic&I  responses  for  the  abc,  abd,  and  at  c capon- 
ants  are  shown  in  vector  diagraa  fora  froa  which  the  total  response 
vector  is  constructed.  Typical  response  for  the  case  where  w  > w 

8 

le  shown  in  the  sketch  below  for  control  gain,  K0,  positive  and 

K  negative. 

8 


The  figure  shows  that  addition  of  coapensstion 
produces  a  system  response  whose  gain  and  phase  are  significantly 
different  froa  the  uncompensated  response.  Regardless  of  whether 
KQ  is  positive  or  negative,  the  maximum  gain  (vector  length)  in¬ 
creases.  For  Kc  positive,  the  originally  stable  system  (vector 

abe)  is  made  unstable  because  of  increased  gain  and  additional 

positive  phase  of  &  .  For  K  negative,  the  marginally  stable 

c 

system  is  made  more  stable  despite  increased  gain  because  of  the 


negative  phase  shift.  The  stabilizing  effect  of  negative  Kq 

agrees  with  the  intuitive  conclusion  that  a  decrease  in  turbine 
speed  in  required  to  offset  the  effects  of  an  increase  in  flow, 
that  acconpanies  increasing  suction  pressure,  P  • 

®  a 

(k)  Structural  Frequency  Less  than  Suction  Line  Fre¬ 


quency.  The  stability  diagram  for  the  case  where  the  structural 
frequency  is  less  than  the  suction  line  frequency  differs  con¬ 
siderably  from  the  previous  case. 


gain  for  Kq  is  again  destabilizing  because  of  gain  increase  and 
unfavorable  phase  shift.  -native  decreases  gain  slightly  and 
shifts  the  phase  of  the  *ponse  into  a  aore  stable  condition. 

It  is  interesting  to  note  that  the  ease  where 
atructuai  frequency  and  suotion  line  frequency  are  equal  corres¬ 
ponds  to  conditions  for  which  the  se  vector  and  the  abd  vector 
ere  equal  in  Magnitude  but  opposite  in  sign  and  hsnes  oanesl. 

The  response  diagram  then  reduces  to  ths  sbe  vector  which  is  ex¬ 
actly  equal  to  the  uncompensated  case.  It  is  also  Important  ts 
note  that  the  magnitude  of  the  abc  response  is  greatest  and  ths 
phase  most  unfavorable  for  the  condition  where  structure  and 
suction  lins  frequency  are  equal.  Thu  under  conditions  which 
lead  to  the  least  stable  configuration,  the  compensation  doss 
little  or  nothing  to  improve  stability. 

d.  Control  with  Suction  Pressure  Feedback.  If  suction 
pressure  instead  of  acceleration  is  used  as  the  control  variable, 


45 


stability  can  bs  improved  in  some  eases  by  an  effect  that  ;  equiv¬ 
alent  to  separating  the  suction  line  frequency  froa  the  stx^tural 
frequency.  With  P#  feedback  through  the  gas  generator/turbine  loop, 

the  P transfer  function  assunes  the  fora  shown  below! 


and  P  is  govern  by 

«m2» 
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(28) 


He  note  here  again  that  when  ■  0  (when  structural 

frequency  and  suction  line  frequency  are  equal)  that  the  K  +  K  /I 

p  tr 

+  ’’B  terms  in  numerator  and  denominator  cancel  indicating  that  the 

compensation  resulting  from  the  control  function  has  been  nullified. 

Thus,  as  in  the  case  of  acceleration  feedback,  the  least  stable 

ease  is  not  isprored  by  presence  of  the  control  compensation. 

If  the  control  gain,  K  ,  is  aade  such  larger  than  K  ,  and  if  r  » 

c  p  s 

1,  which  it  uaually  is  for  the  frequency  range  of  interest,  the 

system  function  can  be  approximated  by  - 
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The  resonant  frequency  of  this  function  is  given  by  the 


equatiomi 


(30) 


The  resonant  frequency  is  lowered  if  K  is  negative  sad 

e 

i 

K  I 

In.  o  .  Thus,  theoretically,  at  least,  it  aopears  possible  to 

v= 

shift  the  resonant  frequency  of  the  propulsion  systes  response 
function  below  any  structural  frequency  that  nay  be  troublesome. 


In  the  numerator  of  the  response  function,  the  effect 
of  the  punp  gain  is  alaost  completely  nullified  by  the  control 
function.  The  total  systen  aaxinua  response  including  structure 
is  shifted  to  a  relatively  stable  phase  condition  as  shown  in 
the  figure  below; 


In  order  to  test  the  feasibility  of  setting  K,  to  nul¬ 
lify  the  punp  gain  and  shift  the  effective  feed  line  frequency, 
the  percent  change  in  gas  generator  flow  required  nust  be  com¬ 
puted.  Paraneter  values  corresponding  to  a  Titan  III  engine 
systea  were  put  in  the  aodel  to  determine  that  about  %  modulation 
of  fuel  flow  to  the  gas  generator  was  required  to  reduce  a  typical 
auction  line  frequency  by  a  factor  of  4. 


4? 


The  large  gain  magnifies  the  effects  of  uncertainties  resulting 
trm  parameter  tolerances.  Thus,  if  this  method  of  control  were 
«Mi,  it  would  he  necessary  to  conduct  extensive  analysis  to  de¬ 
termine  whether  adequate  phase  aargin  could  be  held  for  realistic 
tolerances  placed  on  paraaeters  that  strongly  affect  system  phase. 

2.  SPECIFIC  APPLICATION 

a.  Introduction.  The  application  of  turbine  speed  control 
to  eliminate  FOQQ  is  based  on  the  use  of  the  propellant  puap  as 
an  amplifier  in  which  a  180  degree  phase  lag  is  aaintained  between 
the  discharge  line  flow  and  the  suction  pressure.  The  net  result 
ie  a  saooth  propellant  flow  to  the  main  engines  and  the  associated 
absence  of  any  thrust  perturbation* 

A  simulation  of  pump  turbine  speed  modulation  which  was 
conducted  on  the  analog  computer  wae  based  on  a  model  of  the  Titan 
II  Stage  1  bo  witer  engine  configuration.  The  basic  system  for  thia 
configuration  is  shown  in  block  diagram  form  in  Figure  25.  The  oxi- 
diser  and  fuel .systems  as  shown  begin  with  the  auction  line  followed 
by  the  puap,  discharge  line,  and  injectors  which  supply  propellants 
te  the  combustion  chamber.  The  pumps  are  driven  mechanically  by 
the  turbine  through  reduction  gears  with  energy  for  the  turbine 
supplied  by  a  fuel-rich  gas  generator.  The  propellants  for  the  gas 
generators  are  tapped  from  the  puap  discharge  lines,  fed  through 
the  footstrap  lines  and  metered  through  a  fixe  I  area  cavitating 
venturi.  Thrust  ia  controlled  primarily  by  proper  sizing  of  the 
cavitating  venturi.  The  nominal  fuel  flow  rate  to  the  gas  generator 
ia  about  11.0  lbs/sec  with  an  oxidizer  flow  rate  of  approximately 
1.0  lba/aec.  Since  a  large  part  of  the  propellant  flow  to  the  gas 
generators  ia  fuel,  thrust  can  be  controlled  by  varying  the  fuel 
flow  to  the  gaa  generator  and  replacing  the  presently  used  fixed 
area  cavitating  venturi  with  a  variable  area  cavitating  venturi 
Uf ad  as  the  control  valve  illustrated  in  Figure  25 . 

Theoretical  considerations  indicate  that  gas  generator 
temperatures  and  mixture  ratios  over  a  +  6?6  thrust  range  are  well 
within  design  limits  for  the  particular  eystem  under  consideration. 

b.  Propulsion  System  Equations,  Pump  Through  Injector.  The 
system  equations  which  were  used  to  develop  the  mathematical  model 
for  the  analog  simulation  were  derived  for  a  system  that  includes 
propellant  pumps,  discharge  lines,  and  min  engine  injectors. 

(See  appendix  Ilia  for  those  equations).  The  complete  model  includes, 
analytical  deacriptioaa  of  the  puap,  turbine,  gas  generators  and 
discharge  as  well  <is  bootstrap  feedlines. 


48 


7igur»  2?  Turbin*  Sp**4  Control  Sehomtic 


c.  Model  Development ,  Linearized  Equations.  A  linearization 
of  the  ays tea  equations  for  programming  on  the  analog  computer  was 
performed  for  two  reasons.  First,  the  simulation  was  based  on  small 
variations  about  nominal  or  steady-state  operating  conditions. 

Thus,  the  original  equation  non-linearities  arising  from 
complex  coordinate  coupling  in  the  bootstrap  and  turbine  system  as 
well  as  flow  non-linearities  were  transformed  by  the  well-known 
differential  method  expressed  in  general  as: 


■  t 

.  dx. 

(32) 

dx. 

•  i 

i  ■  1 

i 

The  linearization  made  it  possible  to  study  stability  of 
the  system  by  conventional  frequency  response  methods.  The  method 
used  in  linearizing  the  funotioma  describing  the  basic  propulsion 
system  can  best  be  shown  with  an  example  using  the  oxidizer  pump 
equation  (89)*  The  effect  of  small  changes  of  Pqb,  Nq,  and  y 

cn  (P  .  -  P  )  can  be  functionally  written  as: 
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or,  (P  .  -  P  )  -  K  P  +  K  _  K  +  Km_  ¥ 
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Where:  K 
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All  variables  in  eq,  (  34  )  now  represent  small  changes 
about  a  nominal  set  of  conditions.  These  conditions  or  steady- 
state  values  are  used  to  determine  the  constants  or  coefficients 
in  the  perturbation  equations.  This  procedure  was  used  to  derive 
the  set  of  linearized  equations  which  appear  in  Appendix  Illb. 
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d.  Analog  Modal  of  Complete  Simulation.  The  combined  aqua* 
tiona  of  appendixes  Illb  and  II Ic  comprised  the  complete  system  for  the 
analog  simulation* 

Interdependence  of  these  equations  as  combined  is  demon¬ 
strated  by  the  block  diagram  of  Figure  26.  The  propulsion  sub-systems 
are  outlined  in  the  flow  diagram  according  to  their  basic  functions. 

The  remainder  of  the  system  including  the  structure  and  pump  dynamics 
couple  with  the  propulsion  system  and  are  given  in  detail  in  appendix 
II  Ic. 


(1)  Gas  Generators,  Turbine  Dynamics*  The  response  of 
the  gas  generator,  turbine,  and  pump  pressure  and  flow  to  variations 
in  bootstrap  system  is  dominated  by  a  first-order  time  lag  associ¬ 
ated  with  the  turbine  and  pump  rotary  inertia.  The  lag  in  this  part 
of  the  system  accounts  for  approximately  80*  phase  lag  at  nP090" 
frequencies  of  about  11.0  cps.  This,  of  course,  assumes  no  time 
lag  in  the  gas  generator  combustion  which  although  small  remains 

as  one  of  the  uncertainties  in  the  analysis. 

(2)  Fuel  and  Oxidizer  Pump  Dynamics.  Both  the  fuel  and 
oxidizer  pump  are  represented  by  gains  and  resistances  as  shown  in 
the  equations  and  the  blocks  of  Figure  26.  This  representation  is 
consistent  with  flight  data  which  shows  discharge  and  suotion  pres¬ 
sures  that  are  consistant  with  the  linear  representation. 

(3)  Oxidizer  Discharge  Dynamics.  The  oxidizer  discharge 
system  is  characterized  by  short  discharge  lines  from  the  pump  to 
the  injectors.  The  bootstrap  line  for  the  gas  generator  branches 
off  from  the  discharge  line  and  presents  a  very  high  impedance  to 
pressure  oscillations  is  the  discharge  system.  This  Impedance  is 
in  the  fora  of  a  large  resistance  caused  by  the  relatively  small 
diameter  bootstrap  line  and  cavitating  venturi  represented  by 
(Kg??)  in  Figure  2$.  Consequently,  any  pressure  disturbance  in  the 

oxidizer  discharge  line  is  greatly  attenuated  before  reaching  the 
gas  generator  injector.  The  high  impedance  of  the  gas  generator 
bootstrap  lines  relative  to  the  sain  pump  discharge  lines  make 
it  possible  to  use  the  gas  generator  loop  for  control  without  fear 
of  causing  instability  because  of  the  Inherently  large  tine  lag. 

(4)  Fuel  Discharge  Dynamic e.  The  main  fuel  pump  discharge 
system  contains  regenerative  cooling  tubes  whose  inertance  produces 

a  significant  time  lag  at  POGO  frequencies.  The  presence  of  the 
time  lag  is  first  evident  from  a  comparison  of  coefficients  G„  and 

G^0  (See  appendix  Illb)  in  the  systes  equations.  The  remainder  of 

this  fuel* system  is  much  like  the  oxidizer  system  since  the  fuel 
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bootstrap  line  presents  a  high  impedance  to  the  fuel  pump  discharge 
pressure* 

(5)  Fuel  Controller  and  Qas  Generator.  System  control 
as  previously  stated  was  accomplished  by  replacing  the  cavitating 
venturi  in  the  fuel  footstrap  line  with  a  servo  operated  control 
valve.  The  resulting  analytical  expression  which  describes  the 
behavior  of  the  controlled  fuel  flow  to  the  gas  generator  is  then 
written  as: 

P.  *  K_„  W..  +  K 

fc  023  fb  c 

where  the  primary  feedback  for  control  is  acceleration  at  the 

engine  gimbal  point.  The  gain  of  the  system,  K  ,  and  the  leads 

c 

or  lags  and  r ^  respectively  are  design  parameters  which  are 

optimized  for  the  desired  control.  With  Pfc  remaining  relatively 

constant  due  to  high  impedance  of  the  bootstrap  line,  the  flow 
rate  to  the  gas  generator,  is  essentially  proportioned  to  the 
accelerations,  g^in  with  a  lead  or  lag  as  dictated  by  or  r^* 

Use  of  acceleration  as  the  feedback  variable  rather  than  another 
parameter  such  as  suction  pressure  is  preferred  because  control 
of  acceleration  Itself  is  the  final  requirement. 

e.  Results  of  Simulation.  The  bootstrap  control  system  was 
studied  by  an  open  loop  stability  analysis  using  an  analog  computer. 
With  the  loop  opened  at  acceleration,  g,  frequency  responses  were 
computed  and  compared  with  an  uncompensated  configuration.  The 
uncompensated  or  no-fix  case  shows  marginal  stability  of  about  20° 
phase  margin ^Figure  2?. 

When  the  compensation  was  introduced  by  changing  the  con¬ 
trol  gain,  K  ,  greater  stability  is  achieved.  Although  the  absolute 
c 

value  of  the  gain  at  resonance  has  increased,  the  additional  lag  in 

the  system  increases  the  phase  margin  to  a  value  of  about  98*  for 

the  case  of  It  »  160# 

0 

It  becomes  apparent  that  as  is  increased  beyond  a  value 
of  160,  the  Nyquist  swings  back  into  a  region  of  less  stability  at 
the  higher  frequencies.  The  particular  numerical  value  of  ■  160 

is  not  in  itself  very  significant.  However,  for  the  example  chosen, 
it  shows  that  control  gain  cannot  be  increased  indefinitely  without 
adversely  affecting  stability.  The  effect  of  introducing  a  lead 
or  lag  by  way  of  and  r 2  is  shown  by  part  (d)  of  Figure  27  in 
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which  a  lead  of  about  30*  at  13  cps  has  been  introduced.  A  cor¬ 
responding  loss  of  phase  margin  is  the  result.  Likewise,  if  a 
lag  is  applied  to  the  system,  a  shift  toward  a  more  stable  region 
will  occur  provided  the  gain  is  not  so  large  that  stability  of 
higher  system  modes  is  diminished. 

f.  Other  Analysis  Considerations. 

(1)  Oxidizer  vs  Fuel  Control.  Two  modes  of  control  are 
possible  using  modulation  of  either  fuel  or  oxidizer  bootstrap  flows. 
Fuel  control  as  studied  in  the  preceding  analysis  can  be  replaced 
by  control  of  the  oxidizer  flow.  However,  since  the  steady-state 
oxidizer  flow  rate  is  approximately  10$  of  the  fuel  flow  rate,  the 
oxidizer  control  gains  would  be  much  more  critical.  From  a  prac¬ 
tical  viewpoint  the  parameter  tolerances  on  the  oxidizer  valve 
system  would  require  closer  control  assuming  the  same  percent  change 
in  main  engine  thrust  as  can  be  achieved  indie  fuel  control  system. 

(?)  Gas  Generator  Mixture  Ratio  Insensitivity.  The 
analysis  revealed  that  the  system  is  dynamically  quite  insensitive 
to  gas  generator  mixtures  ratio  because  mixture  ratio  effects  are 
overshadowed  by  the  sensitivity  of  the  turbine  spouting  velocity, 

Vi,  to  changes  in  WQb  and  Sf^.  This  can  be  seen  by  a  comparison 

of  K*  and  with  Kq1  and  Kq2  in  the  equation  of  appendix  Illb. 

However,  this  assumes  that  the  mixture  ratio  of  «ust  remain 


between  fairly  narrow  limits  on  the  very  fuel-rich  mixture  ratio, 
as  seen  in  Figure  86  of  appendix  IHa.  If  the  mixture  ratio  falls 
below  a  value  of  0.06?,  the  C*  falls  off  rapidly.  The  total  effect 
is  t,o  make  turbine  inlet  pressure  very  sensitive  in  this  region  which 
would  cause  sudden  changes  or  non-linear  response  of  the  gas  generator 
combustion  process.  The  P0G0  simulation  placed  a  lower  limit  on 
W  of  .067  so  that  this  non-linear  region  ne.-  not  be  considered. 

An  upper  limit  on  the  mixture  ratio  was  set  at  limit  temperatures, 
below  values  at  which  the  gas  generator  would  burn  up. 

(3)  Mechanical  replacement  of  the  venturi  in  the  gas 
generators  feedline  with  a  control  valve  could  be  done  with  little 
difficulty.  There  appears  to  be  ample  access  to  the  engine  as¬ 
sembly  on  most  liquid  fuel  engines  operating  with  a  bootstrap 
system.  The  valve  to  be  considered  might  be  of  a  variable  venturi 
t.'  pe  but  would  need  lean  and  rich  stops  to  maintain  the  proper  mix¬ 
ture  ratio  range  for  stable  gas  generator  performance. 
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SECTION  V 


FLUIDIC  AMPLIFIER  COMPENSATION 


1.  INTRODUCTION 

The  use  of  a  fluidic  device  or  amplifier  for  control  of 
POGO  la  similar  in  application  to  a  valve  modulating  flow  in 
the  pr  opulsion  system* 

A  fluidic  device  was  studied  because  it  seemed  possible  to 
modulate  flow  without  having  to  drive  sane  mechanical  component 
with  a  sensed  variable  such  as  suction  pressure  or  acceleration. 

Hie  attractive  features  of  a  fluidic  device  in  this  capac¬ 
ity  is  the  absence  of  moving  parts.  The  reliability,  therefore, 
should  be  greater. 


2.  PROPORTIONAL  FLUID  AMPLIFIER 


The  operation  of  a  proportional  fluid  amplifier  can  best  be 
described  with  the  aid  of  the  schematic  diagram  shown  below 


Figure  28  Fluid  Amplifier  Schematic 


Output 

Apertures 


The 
pressure 
Qp.  The 


input  power  to  the  amplifier  is  applied  through  a 
to  the  power  jet,  P  ,  which  in  turn  results  in  a  flow 
control  jets  flows, and  Q  2  combine  with  the  power 
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jet  flow  in  the  interaction  region  to  effect  a  momentum  exchange 
which  deflects  the  output  flow  Q  through  an  angle  0.  The  angle 
of  deflection  is  proportional  to°the  ratio  of  the  mooentums  of 
the  control  stream  and  the  power  stream  as  given  in  the  formula. 


tan  0 


Where:  0 


fluid  velocity 


A  a  nozzle  area 


After  stream  interaction,  the  velocity  profile  of  the  out¬ 
put  stream  is  approximately  gaussian.  The  profile  broadens  and 
decreases  in  centerline  velocity  as  it  moves  downstream  from  the 
interaction  region.  The  flow  splitter  divides  the  stream  which 
is  then  collected  in  the  two  output  apertures.  For  small  angles 
of  deflection,  the  changes  in  output  flows  are  approximately 
proportional  to  the  changes  in  control  flows.  Essential  to  the 
operation  of  the  proportional  amplifier  is  an  interaction 
chamber  as  ihown  which  prevents  the  fluid  stream  from  attaching 
itself  to  the  wall  in  the  output  region  causing  bistable  opera¬ 
tion  of  the  amplifier.  The  flows  gains,  AQ  for  this  type 

°  « 

AQc 

of  amplifier  may  be  as  high  as  10.0  while  pressure  gains  may 
retch  16.0.  However,  the  amplifiers  are  extremely  sensitive  to 
loading  which  can  drastically  reduce  their  gain  as  well  as 
affect  their  operating  stability. 

3.  METHODS  OF  CONTROL  WITH  FLUID  AMPLIFIER 

a.  Fluid  Amplifier/Variable  Resistance  -  Initial  considera¬ 
tions  involved  the  use  of  a  fluid  amplifier  with  its  output  legs 
attached  to  a  converging  Y  joint  as  seen  in  the  schematic 
below: 


^cl 


<*c2 

figure  29  Amplifier/Re si stance 


Basically  this  system  acts  as  a  variable  resistance  when 
placed  in  series  with  another  system  such  as  a  propellant  line. 
The  amplifier  can  route  flow  to  either  leg,  as  shown*  As  more 
of  the  flow  is  routed  through  leg  A,  a  greater  amount  of  fluid 
interaction  occurs  at  the  output  function  which  results  in  a 
greater  pressure  drop  at  this  point.  A  resistance  coefficient 
for  this  function  is  defined  as  follows: 


2  2 


Where:  P  »  fluid  velocity 

Q  *  fluid  flow  rate 


This  assumes  an  angle  of  90*  for  leg  A  with  respect  to 
leg  B  and  the  output  Q^.  The  pressure  drop  is  then  defined  as 


2 

AP  a  (P  P, 

2 


(39) 


for  the  output  function  with  respect  to  flow  past  the  junction* 
Theoretical  values  of  (  vary  with  flow  in  A  and  range  from  0.0 

to  2.0  with  experimental  observations  indicating  the  maximum 
value  to  be  somewhat  nearer  1.7  to  1.8, 
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1)  Series  Suction  Line  Configuration  Method 

a)  One  method  of  using  the  variable  resistance  con¬ 
figuration,  inserts  the  amplifier  in  the  suction  line  as  •  series 
component*  The  accompanying  schematic  of  the  system  shows  the 
major  components. 


This  configuration  shows  one  of  the  control  jets  operating  as 
a  bias  so  that  the  main  stream  flow  is  not  deflected  at  steady 
state  pressures.  As  suction  pressure,  P  ,  begins  to  increase, 

the  output  flow  swings  into  leg  #2.  Since  the  resistance  of  leg 
#2  is  higher  than  leg  A,  the  flow  rate  decreases  thereby  tending 
to  counteract  the  effect  of  increased  pressure.  Ideally  the 
fluidic  amplifier  can  maintain  pressure  at  the  pump  inlet  at  a 
constant  level  for  small  disturbances  in  the  suction  line. 

b)  Practical  Considerations  and  Results*  Two  basic 
problems  arise  with  the  fluidic  amplifier  used  as  a  variable  re¬ 
sistance  configuration.  First,  operation  of  the  power  jet  and 
control  jet  at  the  same  pressure  is  not  practical  or  in  most 
cases  even  possible.  An  alternative  method  controls  flow  in 
the  discharge  line  where  the  operating  pressures  are  suffici¬ 
ently  high.  The  discharge  line  of  a  typical  system  has  pressures 
of  1200  +  20  paia  (steady  state  +  oscillatory)  as  compared  to 
suction  line  pressures  of  45  +  20  psia  (steady  state  +  oscilla¬ 
tory).  However,  the  control  jet  generates  enough  pressure  noise 
to  obscure  the  signal  present  in  the  discharge  line.  The  second 


reason  for  discarding  the  series  suction  line  configuration  stems 
from  the  fact  that  a  20St  pressure  drop  developes  across  the  fluid 
amplifier  even  under  the  aost  ideal  conditions*  A  steady  state 
pressure  drop  of  this  magnitude  degrades  performance  to  such  an 
extent  +  at  it  cannot  be  tolerated  in  the  suction  line  of  pre¬ 
sent-da^  propulsion  systems.  Therefore,  this  particular  method 
of  compensation  was  abandoned. 

2)  Modified  Amplifier  -  Suction  Line  Compensation 

a)  Method.  This  method  was  conceived  as  a  possible 
way  of  introducing  a  variable  resistance  into  the  suction  line 
without  adding  a  steady-state  pressure  drop.  A  diagr^a  of  basic 
elements  of  this  system  follows: 


The  actuator  shown  has  replaced  the  true  fluid  amplifier  although 
its  action  produces  similar  results  by  means  of  the  flow-flipper 
vane  used  to  deflect  the  fluid  flow  into  either  of  two  paths. 

The  high-resistance  path  again  creates  a  virtual  resistance  to 
fluid  flowing  in  the  suction  line  while  the  other  path  is  ef¬ 
fectively  a  j  >ro  resistance  in  the  suction  line.  The  actuator 
as  shown  has  as  its  fluid  power  source  the  steady-state  pressure 
from  the  discharge  system,  P,.  The  maximum  tolerable  flow  that 
can  be  fed  back  without  significant  loss  of  performance  is  about 
556  of  steady  state  -  for  a  typical  system  this  is  about  27.0 
lbs/sec.  As  suction  pressure  increases,  *he  spring  loaded  piston 
drives  a  flipper  vane  forward  resulting  in  more  fluid  flowing  into 
the  high-resistance  leg  of  the  compensator  discharge  system. 
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The  higher  resistance  that  the  suction  iue  now  sees  because  of 
fluid  interaction,  under  ideal  conditions,  will  compensate  for 
the  increased  pressure  in  the  suction  line  to  give  zero  net  flow 
past  the  interaction  region. 

b)  Results.  If  all  of  the  flow  from  the  compensator 
is  directed  through  the  low  resistance  discharge  path  at  steady- 
state  suction  pressure,  a  half  wave  rectified  flow  variation  ac¬ 
companies  the  positive  pressure  pulses.  Typical  values  from  an 
analysis  in  which  idealized  conditions  were  assumed  and  the  fol¬ 
lowing  system  parameters  used,  Indicate  a  compensation  in  the 
amount  indicated  in  the  following  figure: 

P„  ■  45  +  20  paia 

*  900.0  psia  (steady  state) 

<  «  2.0 

2 

A(  *  36.3  in  (suction  line  area) 

Ad  »  0.22  in2 

Q_  -  540.0  Ibs/sec 

S 

Qd  =  27.0  lb*/ see 
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This  indicates  s  sextans  pressure  compensation  of  about  10  psia 
on  the  positive  half  cycle  or  a  reduction  in  the  POGO  content 
of  the  suction  pressure  response  of  about  25%.  Any  compensation 
by  this  method  which  would  work  on  the  negative  half  of  the  cycle 
would  bo  implemented  by  a  shift  in  biasing  or  operating  point 
of  the  compensator.  This  would  also  necessitate  an  additional 
steady  state  suction  line  resistance  which  as  before  ie  intoler¬ 
able.  In  either  case,  the  total  pressure  compensation  cannot 
exceed  10  psia. 

e)  Fluid  Amplifier  Compensation  and  the  Discharge 
System.  As  was  previously  mentioned,  it  is  necessary  for  pres¬ 
sures  to  the  control  jets  of  a  fluid  amplifier  to  be  greater  than 
the  pressure  to  the  power  jet  for  proper  performance.  This  pre¬ 
cludes  the  use  of  the  discharge  system  pressures  for  any  control 

purposes  of  an  amplifier  operating  on  flow  in  a  discharge  system 
As  we  have  seen,  the  signal  would  indeed  be  small  when  compared 
to  the  noise  generated  at  the  control  jet  orifices  end,  therefore, 
made  such  a  system  ineffective. 

4*  CONCLUSIONS  AND  RECOMMENDATIONS 

Fluid  amplifiers  although  attractive  aa  simple  control 
devices  are  generally  impractical  for  POGO  stability  compensa¬ 
tion  because  of  the  wide  range  of  operating  parameters  which 
exist  in  s  modern  rocket  propulsion  system.  Noise  levels,  inter¬ 
nal  stability  and  effects  of  loading  on  performance  could  be 
defined  only  after  full-scale  tests  under  intended  operating  con¬ 
ditions.  If  further  investigations  are  carried  out,  it  would 
be  advisable  to  consider  testing  to  determine  performance  charac¬ 
teristics  with  control  jet  pressures  of  about  the  same  order  of 
magnitude  as  those  of  the  power  jet. 
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SECTION  VI 

CAVITATION  COMPLIANCE  MODEL  DEVELOPMENT 


1.  INTRODUCTION  AND  MODEL  DESCRIPTION 

In  the  analysis  of  structure  -  propulsion  system  coupled 
instabilities ,  turbopump  cavitation  compliance  has  long  been 
subject  to  uncertainty  that  has  limited  the  accuracy  of  stabi¬ 
lity  prediction.  In  order  to  obtain  a  better  understanding 
of  how  the  cavitation  process  affects  the  cavitation  compliance, 
a  mathematical  model  development  study  was  initiated  to  provide 
a  description  of  turbopump  operation.  This  effort  has  been 
carried  over  into  the  present  program  and  is  continuing  at  this 
time.  A  summary  description  of  the  development  and  results 
which  have  seen  achieved  to  date  is  contained  in  the  following 
sections. 


a.  Model  Requirements.  The  objective  of  the  cavitation 
compliance  analysis  was  to  derive  mathematical  descriptions 
that  c  >uld  be  related  directly  to  the  thermodynamic  and  fluid- 
mechanic  conditions  in  a  turbopump.  Semiempirical  approaches 
were  discarded  because  of  their  demonstrated  deficiencies  when 
dealing  with  the  problem  of  pump  head  rise  under  conditions  of 
cavitation. 

The  following  sequence  of  steps,  therefore,  were  con¬ 
sidered  to  constitute  a  logical  analysis  approach: 

1)  Development  of  the  basic  turbopump  flow  equations 
into  which  two-phase  flow  phenomena  could  be  in¬ 
corporated  and  that  could  later  be  expanded  to 
include  more  complex  flow  situations. 

2)  Development  of  a  two  phase  thermodynamic  model 
that  was  independent  of  time  and  conditions  of 
nueleation  and  that  could  be  combined  with  the 
flow  equations  to  give  a  description  of  turbopump 
cavitation  compliance. 

3)  Development  of  a  two  phase  thermodynamic  model 
that  included  time  effects  and  dependence  on  nu.  lei 
and  that,  when  combined  with  the  flow  equations, 
would  also  yield  an  improved  description  of  cavi- 

. ation  compliance. 
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b.  Turboprop  Flow  Equations,  General  Description.  The 
development  of  the  turboprop  flow  equations  began  with  the  as¬ 
sumptions  that  1)  tip  clearance  flow  was  not  present  2)  There  were 
no  vapor  filled  separation  cavities  within  the  impeller,  and 
viscous  effects  were  negligible.  Furthermore,  the  three-dimen¬ 
sional  flow  equations  were  found  to  be  too  complex  for  solutions 
consistent  with  the  scope  of  the  program;  and  in  their  place, 
a  quasi-three-dimensional  potential  flow  approach  was  taken. 

Many  of  the  techniques  used  in  the  development  were  modifications 
of  similar  techniques  used  for  analysis  of  centrifugal  compres¬ 
sors  (Ref.  1  ). 

The  quaai-thre e-dime ncional  analysis  proceeds  first 
with  the  solution  of  the  flow  problem  in  the  meridional  plane. 

It  is  assumed  that  channel  flow  exists  between  the  impeller  blades 
and,  therefore,  a  stream  surface  the  same  shape  as  the  blade  plane 
is  located  midway  between  blades  (Fig.  35  ).  The  relationship 
between  the  stream  surface  and  the  meridional  plane  is  illustrated 
in  Figure  5fc.  The  flow  equations  sure  written  for  the  stream  sur¬ 
face  and  then  referenced  to  the  meridional  plane.  Solution  of 
the  equations  yields  a  description  of  the  flow  field  in  the  meri¬ 
dional  plane  as  shown  in  Figure  3b  .  A  detailed  derivation  of 
the  meridional  plane  flow  equations  is  given  in  Appendix  1  of 
reference  2  • 

The  second  phase  of  solution  involves  rotating  ane  of 
the  meridional  plane  streamlines  and  its  associated  streamtube 
around  the  impeller  axis  to  form  another  stream  surface  and  a 
streamtube  of  revolution.  The  flow  equations  are  written  for 
the  stream  tube  of  revolution,  and  their  solution  yields  a  blade- 
to-blade  description  of  the  flow. 

One  further  simplification  of  the  flow  equations  for 
the  streamtube  of  revolution  was  derived  that  allowed  additional 
reduction  of  problem  time.  The  simplification  is  accomplished 
by  applying  a  number  of  conformal  transformations  to  the  flow 
equations.  The  advantage  of  performing  the  transformations  is 
realised  when  several  different  flow  conditions  in  a  pump  re¬ 
quire  fluid-dynamic  description.  Once  the  coordinate  trans¬ 
formation  is  performed,  all  prop  operating  conditions  are  analyzed 
in  the  simplified  coordinate  system.  Derivations  of  the  stream 
tube  of  revolution  flow  equations  and  the  accompanying  coordinate 
transformations  are  given  in  reference  1  • 
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Solution  to  the  flow  equations  for  both  the  meridional 
plane  and  the  atreastube  of  revolution  is  accomplished  by  numer¬ 
ical  procedures.  Due  tc  the  potentially  large  number  of  flow 
conditions  that  night  require  examination  for  a  given  pump  and 
the  large  number  of  calculations  required  for  a  single  solution, 
the  problems  were  programed  for  the  digital  computer  using  Fortran 
IV  language.  For  this  investigation,  two  programs  were  written. 

The  first  program  performs  the  coordinate  transformations  mentioned 
above.  The  second  program  combines  the  streamtube  of  revolution 
flow  equations  with  a  thermal  equilibrium  cavitation  model  to 
give  a  blade-toblade  description  of  the  compressible  flow  problem. 

c .  Cavitation  Compliance  Problem  Solution.  Using  the  equa¬ 
tion  and  computer  program  development  described  above  the  turbo¬ 
pump  nonseparateu  flow,  thermal  eauilibrium  cavitation  compliance 
is  evaluated  as  follows: 

The  meridional  plane  flow  equations  are  used  elong  with 
the  impeller  geometry  and  pump  inlet  conditions  to  obtain  &  de¬ 
scription  of  the  flow  streamlines  in  the  meridional  plane.  The 
succeeding  step  of  the  analysis  involves  selection  of  a  particular 
streamline  in  the  meridional  plane  and  defining  as  a  function  of 
radius  the  width  of  its  associated  stream  tube  (distance  between 
streamlines  adjacent  to  the  chosen  streamline)  and  angle  a  the 
streamline  makes  with  the  impeller  axis.  These  quantities  are 
tabulated  along  with  coordinates  of  the  impeller  blades  and  con¬ 
stitute  input  data  for  the  coordinate  transformation  program. 

The  coordinate  transformation  program  transforms  the  input  data 
from  the  R,  0,  Z  coordinates  to  the  E,  F  plane  and  subsequently 
to  the  f , n  plane.  An  example  of  the  transformations  is  shown 
in  Figures  35  thru  38  • 

Figure 35  shows  the  pump  blades  in  the  r,  0  plane.  Fig¬ 
ure  37  shows  the  blades  transformed  to  the  E,  F  plane.  Also- 
shown  *,re  the  f  coordinates.  Figure  38  shows  the  blades  in 
the  f,i)  plar.e  where  it  is  seen  that  they  are  parallel,  flat, 
and  of  zero  thickness. 

.  Following  the  coordinate  transformation,  the  two  phase 
flow  problem  is  solved  in  the  f,ij  plane,  utilizing  the  com- 
nressible  flow  program  described  in  appendix  I  of  reference  2 
The  two  equations  that  yield  the  f  ,tj  plane  solution  are: 
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where  the  quantities  U^,  V  ,  b,  r,  and  sin  a  are  defined  in  the. 
reference  and  are  generated  by  fhe  coordinate  transformation  program. 
The  relationship  between  density  and  enthalpy  in  Equation  41  is 
obtained  from  the  thermodynamic  properties  of  the  working  fluid. 

For  pressures  below  the  saturation  pressure  (corresponding  to  the 
fluid  bulk  temperature)  a  relationship  similar  to  that  shown  in 
Figure  39  can  be  obtained.  For  pressures  above  the  saturation 
pressure,  the  flow  is  considered  incompressible. 

Completion  of  the  solution  requires  transforming  the 
streamlines  and  density  lines  from  the  f ,ij  plane  to  the  E, 

F  plane  ana  finally  to  the  B,  0.  Z  coordinates  where  the  relative 
volumes  of  vapor  and  liquid  can  be  assessed  along  with  pump  dis¬ 
charge  conditions, 

2.  oTATUS  OF  ANALYSIS 

The  pump  modeling  phase  of  the  program  is  currently  involved 
with  debugging  of  the  compressible  flow  program  and  generation  of 
density  pressure  relationships  for  various  pump  fluids. 

Following  the  debugging  operation  it  will  be  necessary 
to  checkout  both  the  coordinate  transformation  program  and  the 
compressible  flow  program.  Part  of  the  phase  I  effort  has  been 
devoted  to  compiling  fluid  properties  from  which  pressure-density 
relations  could  be  formulated.  These  functions  will  be  necessary 
input  to  the  compressible  flow  program.  Figures  39  and  present 
the  variation  of  density  with  pressure  for  water  and  N„  0, 
respectively. 
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Both  functional  relationships  are  bases  on  iaentropic  processes 
below  the  saturation  curve. 

Additional  work  has  been  done  during  the  phase  I  pro¬ 
gram  to  try  and  locate  a  pump  impeller  geometry  and  sufficient 
cavitation  test  data  which  can  be  used  for  program  checkout. 

This  work  is  continuing  at  present. 

3.  CONCLUSIONS 

In  order  to  obtain  additional  understanding  of  turbopump 
operation  under  cavitation  conditions  and  at  the  same  time  pro¬ 
vide  means  of  assessing  the  quantity  of  vapor  participating  in 
the  cavit'iion  process,  the  development  of  a  two-phase  flow 
turbopump  model  is  essential.  The  model  should  relate  as  closely 
as  possible  the  thermodynamic  end  fluidmechanic  processes  oc¬ 
curring  in  the  pump  so  that  the  more  sensitive  parameters  in  the 
system  and  their  effect  on  the  POGO  Problem  can  be  recognized. 

The  model  described  herein  is  a  start  in  this  direction;  however, 
due  to  the  broad  range  of  flow  conditions  possible,  the  model 
will  apply  only  to  certain  pump  operating  conditions.  Additional 
improvements  will  be  necessary  to  describe  some  of  the  dynamic 
pump  performance  characteristics  and  their  relation  to  the  suc¬ 
tion  system  continuity  functions.  The  amplitude  limiting  analy¬ 
sis  performed  in  this  program  has  adequately  demonstrated  the 
difference  in  dynamic  and  steady  state  pump  operation  along  with 
the  need  for  a  realistic  continuity  relationship  for  the  suction 
system. 

If  the  gas  injection  concept  of  POGO  control  is  pursued 
to  any  depth,  it  will  become  necessary  to  evaluate  the  effects 
of  injected  gas  in  the  suction  line  and  pump  on  reducing  suction 
line  natural  frequency  as  wel'  as  assessing  itB  influence  on 
pump  performance.  Here  again  a  two-phase  flow  model  will  become 
necessary. 

Based  on  the  many  potential  uses  and  needs  for  a  turbo¬ 
pump  cavitation  flow  model,  it  is  recommended  that  the  develop¬ 
ment  described  in  the  preceding  paragraphs  be  continued  in  the 
phase  II  effort  culminating  with  the  mapping  of  a  turbopurap  and 
correlation  of  test  and  analysis  results. 


SECTION  VII 


STATISTICAL  GAIN  MAfiQIN  STUDIES 


1.  DESCRIPTION 

1%e  primary  result  of  this  study  is  the  development  of  a 
method  which  relates  a  mean  gain  margin  to  a  probability  of 
stability  particularly  as  applied  to  the  POGO  model.  In  the 
mathematical  description  of  any  physical  system,  there  are 
parameters  which  are  inaccurately  known,  or  are  known  only 
within  broad  limits.  Parameter  uncertainties  make  it  diffi¬ 
cult  to  determine  system  stability  when  a  "worst  case"  stabili¬ 
ty  analysis  indicates  instability  and  the  system  is  stable 
using  mean  parameter  values.  For  systems  of  this  kind,  a 
qualitative  measure  of  stability  is  needed.  Gain  margin  is 
a  qualitative  measure  of  stability,  but  is  of  little  use  when 
mean  parameter  values  produce  marginally  stable  results.  It 
seems  reasonable  to  assume  a  more  natural  approach  is  to  com¬ 
pute  an  actual  probability  of  stability.  Then  knowing  this 
number,  it  is  easy  to  predict  from  N  flights,  the  number 
expected  to  be  stable.  For  simple  systems  it  is  possible  to 
compute  a  probability  analytically;  however,  for  complex 
systems,  an  analytic  computation  of  a  probability  of  stability 
becomes  impossible.  Hence  the  analytical  approach  must  be 
abandoned  in  favor  of  a  computer  oriented  Monte  Carlo  technique 
if  an  answer  is  to  be  obtained.  The  Monte  Carlo  approach  for 
computing  the  probability  of  stability  is  conceptually  equiva¬ 
lent  to  flying  a  large  number  of  missiles;  tabulating  those 
which  are  stable;  and  computing  the  probability  of  stability, 
i.e.,  the  ratio  of  the  stable  missiles  to  the  total  number 
flown.  Specifically,  each  system  parameter  which  is  identified 
as  being  uncertain  is  assumed  to  be  either  uniformly  or  normally 
distributed  with  mean  equal  to  the  design  or  nominal  value  and 
variance  determined  either  from  statistical  considerations,  f.g. , 
2  i  N  2 

*  =  -  aean)  or  empirically.  Parameters  which  are 

related  functionally  to  statistical  parameters  are  assumed  to 
be  totally  dependent,  i.e.,  once  the  statistical  parameter  has 
been  chosen  from  a  distribution,  its  value  completely  determines 
the  value  of  the  functionally  related  parameter.  During  each 
iteration  of  a  Monte  Carlo  run,  all  statistical  parameters  are 
drawn  from  their  respective  distributions  and  are  inputed  along 


with  the  deterministic  parameters  into  the  scathematical  model. 
Hie  gain  Margin  is  then  computed,  and  the  system  stability  is 
deterained  and  tabulated.  After  a  suitable  number  of  itera¬ 
tions,  (typically  about  200)  the  probability  of  stability, 

■eon  gain  Margin  and  variance  are  computed  and  printed.  In 
addition,  a  histograa  of  the  gain  Margin  is  plotted.  The  gain 
margin  histograa ,  properly  normalized,  is  a  probability  density 
function  of  the  gain  Margin. 

a.  Objection  to  conventional  gain  Margin 

A  possible  objection  to  the  above  approach  is  that 
the  computed  gain  Margin  bistegran  depends  on  where  the  feed¬ 
back  loop  is  opened.  In  aultiloop  feedback  systems,  the  gain 
■argin  does  indeed  depend  on  where  the  ?  oop  is  opened  and  is 
not  a  unique  indication  of  stability.  However,  it  is  true 
that  a  stable  or  unstable  systen  is  accurately  deterained 
regardless  of  where  the  aeasureaent  is  taken,  i.e. ,  where  the 
loop  is  opened.  Therefore,  as  the  feedback  loop  is  opened  at 
various  places,  the  resulting  gain  Margin  histogram  will  be 
different  but  the  resulting  probabilities  of  stability  remain 
constant.  In  other  words,  the  probability  of  stability 
reaains  an  invariant  function  of  where  the  loop  is  opened, 
whereas  the  gain  Margin  histograa  is  not. 

2.  MATHEMATICAL  MODEL 

Previous  analyses  concerning  P0G0  weru  performed  using 
the  systen  equations  in  second  order  fora.  For  instance,  to 
compute  the  gain  Margin,  the  systen  equations  are  written  as 


Ay  *  f  (42) 


where  y  and  f  are  coluan  matrices  and  A  is  am  nxn  Matrix  with 

typioal  elements  of  the  fora  A  s2  +  B.  .s  +  C  .  By  setting 

ij  ij 

s  a  jm,  the  elements  of  A  become  complex  constants.  A  complex 
matrix  inversion  followed  by  a  coapl'.ex  Multiplication  yields 
the  solution 

y  -  A"xf  (43) 
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A  frequency  sweep  for  the  open  loop  transfer  function  P  /P 

cout  cin 

is  terminated  when  the  phase  of  this  expression  approaches  zero. 
The  resulting  number  is  the  zero  phase  gain  (ZPG),  from  which  gain 
margin  follows  directly. 

a.  Gain  Sensitivity 

To  study  the  variation  of  the  gain  margin  with  respect 
to  certain  critical  system  parameters,  it  is  first  necessary 
to  determine  what  system  parameters  are  critical.  One  possible 
approach  is  to  increment  each  parameter  individually  and  com¬ 
pute  the  change  in  gain  margin.  Those  parameters  which  result 
in  a  large  change  in  the  gain  margin  are  the  critical  parameters. 

A  better  and  computationally  more  efficient  approach  is  to 
compute  the  eigenvalue  sensitivities,  i.e. ,  the  change  in  the 
system  eigenvalues  divided  by  the  change  in  a  system  parameter. 
Such  an  approach,  however,  requires  the  eastern  to  be  represented 
in  first  order  differential  form;  that  iu  the  system  repre¬ 
sented  by  Equation  42  oust  be  rewritten  in  closed  loop  form  as 


x  =  Ax 


(44) 


Then  to  compute  the  sensitivities  of  Equation  44  requires  solv¬ 
ing  Equation  45  . 

,bk  . 

<a.W  (45) 

da  <  *±t 


where  ’  system  eigenvalue 

x^  =  right  hand  eigenvector 

v^  s  left  hand  eigenvector 


and  x,  y  represents  the  conventional  scalar  or  inner  product 
of  two  vectors.  The  derivation  of  Equation  45  is  given  in 
Appendix  IV. 
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Table  I  Typical  Eigenvalue  Sensitivities  for  Simplified  Model 


Eigenvalue 

Sensitivity  with 
respect  to 

Eigenvalue  Sensitivity 
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In  theory  it  is  always  possible  to  transform  Equation  into 
Equation  44  at  least  via  a  computer  program.  This,  however, 
is  unacceptable  for  in  the  execution  of  Equation  45  ,  matrix 

A  is  required  to  be  in  analytical  for*  so  that  the  expression 
dA/da  can  be  evaluated.  It  was  possible  to  analytically  trans¬ 
form  fro*  Equation  41  to  Equation  44  for  a  simplified  POGO 
■odel  (see  Appendix  V)  but  not  for  the  more  complicated  *oael. 
The  sensitivity  analysis  was  performed  with  the  model  -escribed 
in  Figure  41  and  the  results  are  described  below.  Three  systea 
parameters  were  identified  as  being  cost  sensitive.  The 
selection  process  was  arrived  at  by  ccaputing  the  eigenvalue 
sensitivities  of  all  system  parameters  appearing  in  the  simpli¬ 
fied  POGO  model  and  choosing  the  three  most  sensitive ,  which 
in  this  case  were  identified  as  the  mode  A  structural  frequency 
fuel  suction  line  frequency  w  ^  and  the  oxyidizer  suction 

line  frequency  The  simplified  POGO  model  is  described  in 

ten  first  order  differential  equation;  hence,  the  sensitivity 
analysis  relates  to  the  ten  associated  eigenvalues.  The 
results  are  shown  in  Table  I.  Only  seven  eigenvalues  are 
shown  in  Table  I  since  three  of  the  eigenvalues  are  complex 
and  the  sensitivity  of  the  two  complex  conjugate  poles  are 
identical.  Of  the  ten  POGO  system  eigenvalues,  two  are  apt 
to  be  unstable  depending  on  the  configuration  being  studied. 

The  roots  located  at  -*-.542  +  363. 1  are  unstable  and  have  sensi¬ 
tivities  of  i.18  -  j2.5,  f.18  -  jO.OlJ,  and  -  0.01  ♦  j2.07 
with  respect  to  the  parameters,  w  ,  u  and  a>  respectively. 

This  ii.dicates  the  magnitude  and  direction  of  the  change  in 

the  eigenvalue  divided  by  the  change  in  a  system  parameter. 

The  most  critical  parameter  then,  for  tne  unstable  eigenvalue, 

is  w,. 

A 

3.  RANDOM  NUMBER  GENERATOR 

The  random  number  generator  used  in  the  Monte  Carlo  gain 
margin  study  is  based  on  a  recurrence  relationship  involving 
integers.  Two  uniform  distributions  are  generated  on  the 
interval  0  to  1.  A  change  of  variables  is  then  made  which 
transforms  these  two  uniform  distributions  into  a  normal 
distribution  with  specified  mean  and  variance.  The  genera¬ 
tion  oT  the  uniform  distributions  is  discussed  in  Appendix  VI. 

4.  RESULTS  OF  ANALYSIS 

The  dynamical  model  used  for  the  Monte  Carlo  statistical 


Figure  4l  Simplified  POGO  Model 


inhere: 


■  Structural  response  characteristic 


s  Suction  liDe  response  characteristic 
G  a  Discharge  line  response  characteristic 


A  a  Suction  line  area  at  pump  inle* 
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gain  margin  studies  was  developed  during  analysis  of  the  Titan 
lx  and  III  POGQ  problems.  A  T-III  MCL  configuration  was 
chosen  to  illustrate  the  statistical  stability  analysis  for 
the  following  reasons; 


1)  all  necessary  parameter  values  w  -e  readily 
available  from  previous  studies; 

2)  the  results  of  the  study  are  useful  for  the 
MOL  program; 

3)  much  of  the  data  taken  from  Titan  II  flights 
is  applicable  for  determining  the  means  and 
variances  for  Titan  III  parameters. 


The  Titan  III  POGO  math  model  is  written  in  the  form  of 
Equation  44  with  the  dimension  of  the  A  matrix,  in  the 
uncompensated  configuration,  equal  to  13.  The  uncertain 
quantities  which  may  be  inputed  into  the  A  matrix  at  each 
iteration  of  a  Monte  Carlo  run  are:  gA,  gg,  gc,  gQ  (mode 

A,  B,  C,  D  structural  gain),  “  g»  At  B,  ci 

D,  structural  frequency),  f  c,  fD,  (mode  A,  B,  C,  B 


damping  ratio),  NPSHO  (net  positi/e  suction  head  of  oxidizer, 
and  NPSHF  (net  positive  suction  head  of  fuel).  The  model 
developed  for  the  present  study  treats  pump  gain  (K  )  and 

suction  line  frequency  (ui  )  as  functions  of  NPSH. 

6 


For  typ_cal  values  of  the  means  and  variances  imputed  into  a 
Monte  Cano  run,  see  Table  u  .  Variables  with  zero  deviations 
ax  ;  set  equal  to  their  mean  value.  Since,  in  most  cases,  flight 
data  is  not  available  for  structural  gains,  frequencies,  and 
damping  ratios,  a  heuristic  approach  was  used  in  determining 
their  atatitistical  properties.  Mean  values  are  taken  as  the 
design  value  and  variances  were  determined  by  using  estimates 
of  prediction  accuracy  for  structural  parameters  such  as 
frequency  ana  gain.  All  atatitistical  variables  used  in  this 
study  were  assumed  to  be  gaussian  or  normally  distributed. 

The  option  is  available  however  to  use  uniform  distributions 
if  desired.  The  mean  and  variance  for  the  NPGH  fuel  and  oxi¬ 
dizer  variables  are  determined  from  the  equations 
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ff 

It  e  “  L  NPSH.  (  o) 

i=l  1 

N 

ff2  =  “rL  (NPSH  -^}2  (4?) 

i=l  1 


and  may  be  updated  as  more  flight  aita  becomes  available. 

a.  Typi  cal  Analysis 

Typical  FOGO  stability  analyses  are  performed  using 
’’best  estimates"  for  input  quantities.  When  there  is  a  large 
uncertainty  for  a  particular  input  quantify  it  is  chosen  in  a 
manner  which  minimizes  the  gain  margin.  This  conservative 
approach  is  ^satisfactory  except  when  trying  to  determine 
the  least  possible  system  stability.  For  example  an  analysis 
which  predicts  stability  using  average  input  values  and  in¬ 
stability  using  worst  case  input  values  leaves  unsettled  tne 
fundament* l  question:  "Will  che  vehicle  oscillate?"  A 
quantitative  judgment  of  the  significance  of  gain  margin  can 
be  made  by  treating  the  input  parameters  as  /anuom  variables 
with  specified  distriDutions,  Then  the  gain  margins  computed 
from  ‘■hose  input  parameters  are  in  turn  random  variables  and 
their  associated  probabilities  are  the  desired  quantitative 
results.  In  most  realistic  problems,  a  closed  form  solution 
is  impossible;  and  Monte  Carlo  technique  must  be  used  to 
obtain  the  probability  of  stability.  The  use  of  a  Monte 
Carlo  technique  assumes  that  an  infinite  population  can  be 
adequately  approximated  by  a  finite  population.  This  then 
gives  a  clue  in  determining  when  the  number  of  samples  are 
sufficient  largc.  Determining  an  adequate  number  of  samples 
is  done  by  doubling  the  number  of  iterations  and  noting 
whether  the  answer  changes  significantly.  After  several 
tests  it  was  determined  that  200  iterations  appeared  to  be 
sufficient  to  obtain  reliable  probability  fstimates. 

b.  Computer  Time 

In  any  iterative  scheme  programmed  on  a  digital 
computer  a  few  seconds  saved  per  iteration  results  in  a 
large  time  saving  per  problem.  Considerable  effort  was  ex¬ 
pended  improving  the  convergence  time  required  to  obtain  the 


zero  phase  gain.  The  result  was  a  reduction  in  convergence  time 
from  2  sec  to  0.2  sec  per  zero  phase  gain  calculation.  This  is 
reflected  in  a  computer  time  saving  from  400  sec  to  40  see  p«r 
Monte  Carlo  run.  Of  course,  as  computers  are  improved  and  the 
arithmetic  processors  become  faster,  the  time  per  computation  is 
also  reduced;  in  fact  five  years  ago  the  same  problem  would 
have  taken  more  than  10  minutes  compared  witt.  40  sec  today. 

c*  Results  of  Typical  Run 

The  main  results  obtained  from  a  typical  statistical 
gain  margin  Monte  Carlo  non  are; 

1)  a  histogram  of  the  zero  phase  g.'dLn;  upon  normali¬ 
zation  the  histogram  becomes  the  zero  phase  gain 
density  function; 

2)  the  probability  of  stability; 

3)  the  sample  mean  of  the  zero  phase  gain; 

4)  the  staple  ''-rience  of  the  zero  phase  gain; 

3/  a  lisciug  of  the  computed  zero  phase  gains. 

A  histogram  is  a  plot  of  the  number  of  the  zero  phase  gains 
contained  in  a  specified  interval.  The  abecias>a  of  the  plot 
is  subdivided  into  N  intervals  and  the  number  of  zero  phare 
gains  contained  in  each  interval  is  tabulated  and  plotted  in 
the  form  of  a  bar  graph.  Figure  42  displays  a  typical  ZT1 
istogra*  for  the  Titan  III  vehicle  in  the  MOL  uncompensate.! 
configuration  at  i  *  80  sec.  The  interesting  aspect  of 
Figure  42  t  which  is  true  of  all  histograms  computed  (over  a 
hundred),  is  that  it  is  severely  skewed  with  a  long  tall  from 
the  left  and  very  little  from  the  right.  Obviously  the  zero 
pnase  gain  density  function  is  not  normal  vuvu  though  dll  the 
input  quantities  are  normal. 

If  the  zero  phase  gain  histogram  were  symmetric,  the  mean  of 
the  zero  phase  gains  and  the  zero  phase  ga.  n  obtained  by  using 
the  mean  input  value  should  coincide.  Because  the  zero  phase 
gain  histogram  is  not  symmetric,  the  mean  zero  phase  gain  and 
the  zero  phase  gain  of  mean  inputs  are  not  identical.  Skew¬ 
ing  to  the  left  will  cause  the  mean  zero  phase  gain  to  be 
less  than  the  zero  phase  gains  of  mean  inputs  and  hence 
reflects  a  more  optimistic  picture  of  vehicle  stability  than 
would  otherwise  be  determined. 
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Tne  mean  zero  phase  gain  and  me  zero  pna.se  gain  of  mean  inputs 
are  compared  for  several  configurations  corresponding  to 
flight  time  of  3C-l40  seconds. 


Table  III  A  Comparison  of  Zero  phase  Gains 


Histograms  corresponding  to  flight  times  90  sec  through  140  seconds 
are  shown  in  Figures  42,  4j5 ,44  ,  45  ,  46  and  47,  These  figures 
reveal  that  every  zero  phase  gain  histogram  is  sxewed  to  the 
left  and  each  succeeding  histogram  is  displaced  farther  to  the 
right,  i.e. ,  the  vehicle  tends  to  be  less  stable  with  increasing 
flight  time.  Tt  is  noticed  that  the  probability  of  instability 
increases  rapidly  from  t  =  80  sec  to  t  =  90  sec  even  though  the 
mean  zero  phase  gain  changed  less  than  0.6  of  a  db.  The  reason 
is  obvio  is  upon  examining  Figures  42  and  45  ;  since  the  density 
function  is  skewed,  a  small  shift  to  the  right  greatly  increases 
the  a:*ea  to  the  left  of  the  0  db  line. 

d.  Effects  of  Shifts  in  Parameter  Means 

The  critical  system  parameters  have  been  identified, 
their  distributions  approximated  and  a  Monte  Carlo  gain  margin 
analysis  performed.  A  natural  question  remaining  concerns 
what  effect  the  shift  in  the  mean  values  of  critical  parameters 
has  on  the  zero  phase  gain  distribution  shape  and  mean  value. 
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Accordingly,  a  Monte  Carlo  analysis  was  run  with  the  mean 
values  of  u  ,  NPSHO  and  NPSHF  incremented.  The  results,  for 
a  change  in  the  mean  NPSHO,  are  shown  in  Figures  49»ad  50  . 
Figure  49  shows  the  zero  phase  gain  histogram  with  the  mean 
NPSHO  10#  low  whereas  in  Figure  50  the  mean  NPSHO  is  10#  high. 
For  a  decrease  in  mean  NPSHO  the  mean  zero  phase  gain  increased, 
i.e. ,  the  system  becomes  more  unstable;  as  the  mean  NPSHO  is 
increased  the  mean  zero  phase  gain  is  decreased  and  the  system 
stability  is  increased.  Also  the  spread  or  variance  of  the 
zero  phase  gain  density  function  increases  with  decreasing 
NPSHO.  A  possible  explanation  is  that  as  the  mean  NPSHO  in¬ 
creases,  the  separation  betwet.  'tructural  and  suction  line 
frequencies  also  increases  whicn  results  in  an  increased 
stability  margin.  Also  the  variability  or  the  second  moment 
of  the  zero  phase  gain  density  function  seems  to  be  highly 
dependent  upon  the  separation  of  the  structural  and  suction 
line  frequencies.  That  is,  as  the  suction  line  and  structural 
frequencies  tend  to  become  equal,  a  small  change  in  a  system 
parameter  results  in  a  large  change  in  the  zero  phase  gain; 
whereas  for  the  case  where  the  structural  and  suction  line 
frequencies  are  widely  separated,  a  small  change  in  a  system 
parameter  results  in  a  small  change  in  the  zero  phase  gain. 

The  significance  of  the  above  observation  relates  to  the 
desirability  of  maintaining  a  separation  between  the  structural 
and  suction  line  frequencies  when  system  predictability  is  con¬ 
cerned. 


e.  Skewing  of  the  ZPG  Density  Function 

By  viewing  Figures  42  through  48  ,  the  one  consistent 
characteristic  evident  in  all  the  curves  is  a  skewing  to  the 
left.  An  explanation  evolves  from  the  functional  relationship 
between  the  inputed  random  variables  and  the  zero  phase  gain. 
Functionally  we  have; 

Zero  Phase  Gain  =  f(w.,  NPSHO,  NPSHF) 

A 

Clearly,  since  a>^,  NPSHO  and  NPSHF  are  random  variables,  the 

zero  phase  gain  is  also  a  random  variable.  The  problem  which 
is  solved  by  the  Monte  Carlo  technique  is  the  determination 
of  the  zero  phase  gain  density  function.  If  the  function*! 
relationship  expressed  by  Equation  '48'  were  a  simple  one, 
it  would  be  possible  to  derive  the  distribution  for  the  zero 
phase  gain  and  the  Monte  Carlo  simulation  would  not  be  necessary. 
This  however  is  not  the  case.  The  question  to  answer  next  ia: 


(48) 
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Vhat  kind  of  functional  relationship,  expressed  by  Equation  48  , 
is  required  to  result  in  skewing  of  the  zero  phase  gain  density 
function?  It  can  be  shown  that  to  obtain  etc  swing,  Equation  48 
definitely  met  not  be  linear. 

It  is  possible,  but  tedious  to  construct,  by  a 
grapaical  toohnique,  the  functional  relationship  described  by 
Equation  48  when  the  fora  of  the  zero  phase  gain, u  A,  NPSBO 

and  IfPSHT  density  functions  are  known. 

f.  Plotting  Histograa  rs  Gain  in  db  vs  Plotting  in  Gain  Units 

There  are  two  ways  of  plotting  the  zero  phase  gain  histo¬ 
graa*  the  gain  can  be  plotted  as  a  siaple  ratio  of  input  over 
output  or  the  gain  can  be  plotted  in  teraa  of  db.  The  result  is 
•iaply  a  napping  froa  one  variable  to  another.  For  instance, 
suppose  the  zero  phase  gain  density  function  were  expressible  as 
an  analytic  curve  of  gain  units  y,  i.e. , 


SPG  »  f(y)  (49) 


and  since  db  ■  20  log  y  the  expression  for  the  zero  phase  g»in 
density  function  in  tens  of  db  is  siaply 

ZPG  .  f(10dV/2°)  .  g(db). 


Hence,  it  is  a  siaple  aatter  to  transforn  froa  the  one  coordi- 
**te  ayatea  to  another.  Figure42  which  is  the  zero  phase  gain 
density  function  plotted  as  a  function  of  db  is  replotted  as 
a  function  of  gain  units  in  Figure  51  •  Soac  of  the  ainor 
details  are  different,  but  the  overall  characteristics  regain 
unchanged.  Obviously  the  above  change  of  variable  does  not 
change  the  probability  of  stability,  but  does  change  both  the 
mean  and  variance. 

g.  Probability  Contours 

Once  the  zero  pftase  gain  density  function  has  been 
obtaiaed  it  is  possible  to  coopute  a  probability  of  stability. 
Then  as  conditions  change  with  flight  tiae,  a  new  zero  phase 
gain  density  function  can  be  coaputed  and  also  a  corresponding 
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probability  of  stability.  Such  a  procedure  can  result  in  a 
plot  of  probability  of  stability  versus  fl’ght  tine;  for 
various  flight  tiaes  there  corresponds  various  probabilities. 
Oace  the  *ero  phase  gain  density  functions  have  been  obtained 
it  is  possible  to  specify  requirements  of  a  compensator  capable 
of  aaintaioing  a  particular  probability  of  stability  over  a 
specified  flight  tiae.  For  instance  consider  a  zero  phase 
gain  density  of  the  fora  shown  below  in  Figure  52  . 


Figure  52  Typical  ZPG  density  Function 

Assume  for  purposes  of  illustration  that  the  probability  of 
stability  is  0.25.  It  Is  desired  to  compensate  this  system  so 
that  the  probability  of  stability  is  increased  to  0.75.  By 
integration  it  is  possible  to  determine  the  distance  X  so  that 
795*  of  the  area  of  the  curve  is  to  the  left  of  the  dotted  line. 
The  value  of  X  then  is  the  amount  of  comp  »sation  required  to 
maintain  the  probability  of  utability  at  the  desired  level. 

For  such  an  approach  to  be  valid  we  have  ioplicitly  made  several 
assumptions.  First,  we  have  assumed  that  the  compensator  will 
uot  change  the  shape  of  the  zero  phase  gain  density  function; 
and  secondly,  we  lave  assumed  that  the  compensator  is  in  series 
with  the  rest  of  the  systems  where  the  loop  is  opened.  Under 
these  assumptions,  the  approach  is  valid. 

Using  ine  above  approach,  and  the  zero  phase  ^ain  curves  shown 
in  Figures  42  through  48  ,  plots  of  constant  probability  were 
obtained.  The  three  curves  shown  in  Figure  53  correspond  to 
probabilities  .95»  .75  and  0.5.  The  ordinate  of  Figure  53 
corresponds  to  the  compensation,  in  db,  required  to  saintain 
the  associated  probabilities.  Because  of  the  skewed  nature 
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of  the  density  functions  nil  three  curves  of  Figure  53  tend 

to  b#  grouped  together. 

5.  A  COMPARISON  OF  MONTE  CARLO  METHOD  WITH  RSS  METHOD 
a.  RSS  background 

When  the  coapananta  of  ary  system  are  operating  in 
cassette,  each  may  be  generating  a  redoes  error.  If  the  causa 
©f  tu;.  error  in  one  component  is  unrelated  to  the  cause  of 
:h«  error  in  some  other  component,  they  are  said  to  be  inde¬ 
pendent.  Thus  while  t  “  error  in  one  component  may  be  positive 
the  error  in  another  component  can  be  negative,  etc.  However, 
when  it  ia  necessary  to  add  up  nil  the  ?rrora  to  determine  the 
total  error,  straight  addition  of  1-sigma  values  would  give  a 
peeaimietic  description  of  the  probable  error.  Thus  the  usual 
technique  of  cofsbinlig  1-sigma  ei  ro rn  vh«n  they  are  inde¬ 
pendent  of  coe  aacther  is  the  r >ot~«ym-square  technique. 

If  w^,  a  2»  •  »  *,<rt  represent  several  errors  combining  to 

cans*  an  overall  error,  then  the  overall  error  is  called  the 
1-cigma  ESS  (r^ot-soa-aquare)  and  is  obtained  by: 


ffRSS  * 


‘*i> 


(o,)2  * 

d 
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In  the  gain  margin  version  of  the  RSS  technique  the  system 
parameters  NP.SHO,  NPSHF,  *  4  and  g  are  sequentially  varied. 
The  difference  between  the  'nominal  gain  margin  and 


In  the  gain  margin  version  of  the  RSS  technique  the  system 
parameters  NPSHO,  NPSHF,  <i  ^  and  g^  are  sequentially  varied. 

TSam  difference  between  the  nominal  gain  margin  and  the 
Incremented  gain  margin  is  computed  as  each  parameter  is 
varied,  end  these  are  combined  according  to  Equation  31 
This  number  then  is  estimate  of  the  zero  phase  gain  standard 
deviation*  It  is  interesting  to  compare  this  estimate  of  the 
aero  phase  gain  standard  deviation  with  that  obtained  by  the 
Monte  Carle  technique,  tfsing  the  RSS  technique  on  the  Titar 
III  POGO  model,  the  computed  value  of  a  was 


0 

RSS 


0.4C4 


(51) 


(52) 


Using  the  same  mathematical  modal  in  a  Monte  Carlo  simulation, 
th$  resultant  value  of  a  was  computed  to  be: 
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0.114 


(53) 


Monte  Carlo  * 


Sine*  the  Moot*  Carlo  technique  is  a  simulation  of  the  statisti¬ 
cal  processes  involved,  it  is  reasonable  to  amsume  that  the 
zero  phase  gain  deviation  (ffHoot#  Carlo ^  80  should  b* 

reasonably  accurate.  Hence  the  ESS  estimate  of  the  standard 
deviation  is  overly  pessimistic  in  its  approximation  to  the 
true  parameter. 


CONCISIONS 


The  primary  result  of  this  study  has  been  the  development 
of  a  computer  program  which  relates  a  mean  gain  sargin  to  a 
probability  of  stability.  Any  linear  constant  coefficient 
feedback  system  in  which  the  parameters  can  be  described  as 
random  variables  with  specified  distributions  is  an. enable  to 
this  approach.  Using  the  Titan  III  F0Q0  Model,  several  illustra¬ 
tive  examples  were  worked.  In  particular  the  probability  of 
stability  was  computed,  a  zero  phase  gain  histogram  constructed, 
and  the  sample  mean  and  sample  variance  oi  the  zero  phase  gain 
calculated.  It  was  shown  that  the  shape  of  the  zero  pha*" 
gain  histogram  is  related  to  where  the  measurement  is  taken 
(i.e. ,  where  the  loop  is  opened)  whereas  the  probability  of 
stability  is  not.  The  zero  phase  gain  histogram  when  properly 
normalized  becomes  a  probability  density  function,  Beery 
histogram,  constructed  from  the  Titan  III  FOOD  Model,  displayed 
a  skewing  to  the  left.  This  was  assumed  to  have  been  the  result 
of  the  functional  relationships  between  the  statistical  para¬ 
meters  and  the  deterministic  parameters.  As  more  flight  data 
is  obtained,  it  is  always  desirable  to  use  this  data  to  improve 
previous  parameter  estimates.  This  approach  is  amenable  with 
the  Monte  Carlo  simulation.  With  increased  flight  data,  the 
statistical  input  parameters  (mean  and  variance)  can  he  updated 
using  Equations  46  giving  improved  estimates  of  system 
stability. 

a.  Disadvantage  of  Method 

A  possible  drawback  of  the  Monte  Carlo  method,  wham 
extremely  small  or  accurate  estimates  of  the  probability  are 
desired,  is  that  many  iteration®  are  required.  If  for  Instance 
a  probability  of  stability  of  0.999  is  required,  the  number  of 
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Monte  Carlo  iterations  to  obtain  this  accuracy  would  be  approxi¬ 
mately  5*000*  With  present  computing  speeds,  5POO  iterations 
results  in  approximately  16  minutes  of  machine  time.  To  obtain 
estimates  of  the  extremes  or  tails  of  distributions  without 
taking  many  Iterations,  it  is  possible  to  use  special  sampling 
techniques  which  accentuate  these  extremes.  Special  sampling 
techniques,  however,  were  not  used  in  this  study.  It  seems  that 
the  association  of  a  probability  to  quantitatively  describe  the 
degree  of  stability  is  both  a  natural  and  useful  approach.  The 
price  one  pays  for  this  approach  is  that: 

1)  An  analytical  computation,  in  most  cases,  is 
impossible. 

2)  A  large  amount  of  computer  time  may  be  required. 

The  last  objection,  however,  will  considerably  diminish  as 
computational  speeds  are  increased. 
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SECTION  VI n 


AMPLITUDE  LIMITING 

1.  INTRODUCTION  -  STATEMENT  OF  PROBUM 

In  trying  to  augment  the  basic  understanding  of  the  POGO 
phenomena,  one  hau  to  contend  with  the  observation  that  the  longi¬ 
tudinal  oscillations  on  Titan  II  VeL^ilea  increased  to  a  Uniting 
value  and  then  diainished  to  a  such  ltver  level  or  disappeared  be¬ 
fore  first-stage  burnout.  In  the  first  few  flights  on  Titan  II , 
the  duration  of  instability  as  shown  in  Figure  54  covered  approximately 
the  sane  tine  in  flight;  however,  the  Halting  anplitude  and  the  ampli¬ 
tude  envelope  did  vary  substantially  over  these  flights  which  exhibited 


Figure  54  Titan  II  Structural  Acceleration 

a.  Observed  Variations .  Those  parameters  which  are  considered 
to  be  input  variables  for  the  aath  model  and  were  observed  to  vary 
from  flight  to  flight  are  listed  in  the  table  below.  The  raoge 
of  variation  of  these  parameters  presents  a  formidable  problem  to 
the  investigator  trying  to  model  the  POGO  phenomena  in  a  manner 
capable  of  duplicating  individual  vehicla  behavior. 


TABLE  IV  INPUT  PAR AMI  ER  VARIATIONS 


At  90  Sec  Into  Flight  1 

Miaaile 

Tan's  Top  Pressure 

Steady  State 
Suction  Pressure 

WEI9M 

mm 

N-l 

20.0 

21.5 

66.0 

27.3 

N-4 

16.7 

23.5 

65.5 

30.0 

N-5 

17.2 

19.0 

70.0 

25.9 

N-6 

17*5 

19.5 

71.0 

26.0 

b.  Objectives.  The  objectives  of  tha  amplitude  limit  study 
wars  to  identify  and  dafina  i at  mathematical  fora  of  ayatea  non- 
liaaaritiaa  reapooaibla  for  the  salient  features  of  tha  atabla  POGO 
liait  cycle*  To  this  and,  flight  data  was  used  in  conduction  with 
an  existing  linear  aodel  to  derive  a  ayatea  shose  critical  ayatea 
eleaenta  ware  determined  by  response  magnitude • 

2.  PROPOSED  APPROACH 

Linear  math  models  of  the  coupled  propulsion  ayatea  and 
structure  have  been  successful  in  predicting  the  period  of  instab¬ 
ility  hut  cannot  duplicate  the  amplitude  limiting  that  has  been 
observed  from  the  Titan  II  flights.  Although  the  depletion  of 
propellant  has  an  appreciable  effect  on  both  the  vehicle's  struc¬ 
tural  behavior  and  the  feed  system  response  (Fig*  91  &  92),  the 
shift  in  natural  frequencies  of  the  various  subsystems  is  not 
sufficient  to  confine  the  oscillations  with  the  linear  aodel  alone* 
For  this  reason,  it  became  necessary  to  as suae  the  observed  ampli¬ 
tude  limiting  was  the  result  of  a  nonlinear  phenomenon  somewhere 
within  the  vehicle* 

a*  Possible  Nonlinearities.  Upon  examining  the  steady-state 
governing  equations  used  in  the  linear  POGO  acd«  ,  it  was  possible 
to  see  a  number  of  nonlinearities  that  could  conceivably  enter  in¬ 
to  the  systems  response  for  large  amplitude  oscillations.  Non- 
linearities  oould  appear  within  the  structural  behavior  of  the 
vehicle  aaking  the  second  order  differential  equation  used  in  the 
linear  description  inadequate  for  large  aaplitude  oscillations. 
Seference  7,  for  example,  points  out  that  on  Titan  II  vehicles 
having  low  ullage  pressures  the  region  of  the  stage  I  oxidizer 
and  fuel  tanks  can  be  seen  to  exhibit  buckling  at  certain  tiaea 
in  the  flight  with  consequent  abrupt  loss  of  stiffness.  The 
inability  of  the  propellants  to  sustain  any  appreciable  tensile 
force  causes  a  number  of  other  nonlinear  behavior  patterns  to 
evolve  when  the  fluid  is  made  to  pulsate  at  sizeable  amplitudes. 

In  addition  to  the  nonlinear  fluid  mechanic  aspects,  large  per¬ 
turbations  in  the  propellant  flow  rates  can  alter  the  mixture 
ratio  and  result  in  nonlinear  changes  in  combustion  ana  tnrust. 

Same  of  the  above  nonlinearities  were  investigated  by 
Goldman  (Ref.  6  )  but  only  to  determine  their  effect  on  open 
loop  aystem  gain.  As  pointed  out  earlier,  the  objective  here 
was  to  determine  possible  non-linearities  which  could  cause  a 
closed  loop  system  that  ia  basically  unstable  to  limit  cycle. 

b.  Dominant  Nonlinearities.  It  was  found  that  the  dominant 
nonlirr  i.r*ti«s  are  those  associated  with  the  propellant  feed  system, 


104 


the  turbopump,  and  its  related  cavitation  region*  An  earlier  study 
(Ref.  2)  investigated  the  iaportance  of  this  cavitation  region  and 
how  the  presence  of  vapor  in  the  feed  eye tea  affected  the  system's 
response  to  perturbations.  It  was  in  this  study  that  the  relation 
of  cavitation  coapliance  to  FOGO  amplitude  limiting  was  discussed 
(Section  IV,  6a). 

3.  NONLINEAR  MATHEMATICAL  MODEL 

The  equations  used  to  describe  the  nonlinear  POGO  response 
are  discussed  in  Appendix  VIII  * 

a*  Nonlinear  Model*  The  linear  sath  nod el  used  in  earlier 
studies  was  revised  by  incorporating  two  nonlinear  relationships* 
One  dealt  with  the  pressure  aaplification  through  the  turbopuap 
by  relating  the  suction  and  discharge  pressures ;  the  other  in¬ 
volved  description  of  flow  continutity  between  the  suction  and 
discharge  flows*  Although  the  two  are  inherently  the  result  of 
the  ease  cavitation  phenoaenon,  they  were  considered  independently* 
During  the  studv.  one  relationship  was  held  constant  while  the 
other  was  varied  in  an  atteapt  to  get  soae  idea  of  their  relative 
iaportance. 

b.  Cavitation  Coapliance.  The  mathematical  models  of  the 
cavitation  coapliance  are  derived  and  discussed  in  detail  in  Ap¬ 
pendix  IX  •  In  brief,  the  cavitation  region  was  represented  three 
different  ways. 


(1)  Amount  of  Vapor  Involved.  Assume  a  fixed  mass  of 
vapor  is  present  in  the  cavitation  region  at  any  tiae  during  the 
pressure  oscillations. 

(a)  Adiabatic  Process.  Further  assume  the  fixed 
mass  of  fluid  undergoes  expansion  and  compr  -aicm  without  the  ad¬ 
dition  or  deletion  of  heat.  With  the  above  two  restrictions  the 

following  continutity  expression  between  and  V.  recalls: 

•o  do 


whert  AV 


■*  «  p  , 

ito  ref 


Vdt 


. ...I... 

1 1*4 

if 


and  in  a  similar  manner  for  the  fuel  system 
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1  '*•  . 
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(55) 


(b)  Isothermal  Process.  The  second  representation 
which  was  examined  assumed  the  expansion-compression  process  took 
placo  while  maintaining  a  constant  temperature. 
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(2)  Amount  of  Vapor-Liquid  Involved.  A  third  model  of 
the  cavitation  region  assumed  a  fixed  mass  of  a  liquid-vapor  mix¬ 
ture  to  be  involved  in  the  pressure-volume  changes.  During  oscil¬ 
lations  in  pressure  within  the  propellant  lines*  the  mixture  was 
allowed  to  change  quality.  For  example,  during  that  part  of  the 
cycle  where  the  pressure  was  increasing,  vapor  was  allowed  to  con¬ 
dense  into  liquid.  The  heat  of  vaporisation  released  during  the 
phase  change  was  stored  within  the  fixed  aass  of  mixture  as  random 
kinetic  energy.  This  process  was  assumed  to  follow  a  constant 
enthalpy  line  under  the  saturation  curve  for  the  preliminary  study 
because  of  the  ready  availability  of  temperature  -  enthalpy  re¬ 
lationship  for  the  oxidizer.  (Ref*  8  ). 


(3)  Representative  Curves.  A  typical  plot  of  each  of 
the  three-proposed  models  is  presented  in  Figure  55  .  Each  of  the 
models  had  two  points  in  common: 


(a)  The  cavitation  region  is  occupied  by  a  fixed 
aass  of  fluid  which  undergoes  this  pressure- 
volume  change. 

(b)  Thermodynamic  equilibrium  exists  throughout 
the  cycle;  thus,  vaporization  and  condensation 
are  considered  as  instantaneous  processes 
without  any  associated  time  lag. 
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Although  different  approaches  were  used  in  arriving 
at  the  three  representations,  there  is  surprisingly  little  differ¬ 
ence  in  the  resulting  P^  vs  AW  relationships  as  shown  in  Figure  55 

e.  Turbo pump  Characteristic  Curve.  The  other  nonlinearity 
which  was  used  in  conduction  with  one  of  the  above  forme  of  the  con¬ 
tinuity  equation  was  the  turbo pump  characteristic  curve.  Test  data 
derived  from  steady-state  operation  of  the  puaps  were  used  to  re¬ 
present  the  dynamic  behavior.  Later  in  the  program  a  more  thorough 
examination  of  the  pump  model  was  undertaken  and  results  show  that 
dynamic  pump  operation  departs  significantly  from  the  static  charac¬ 
teristic. 


(1)  Governing  Equation,  The  discharge  pressure  was 
expressed  as  a  function  of  the  auction  pressure  and  flow  rate 
throagh  the  pump. 


»«  <V  -  *»  »i 


where:  R  *  pump's  internal  resistance 

P 

■  discharge  pressure  from  steady-state  head-rise  data 

A  m  discharge  flow  rate. 

wd 

(2)  Linearised  Pump  Model.  Two  different  pump  models 
evolved  from  Equation  58.  The  first  expressed  as  a  linear 
function  of  the  suction  pressure. 

P.  -  (M  ♦  UP  -  H  W, 

«  spa 


where  (M  ♦  1)  ■  lia««rin*d  pump  gain. 


A  typical  plot  of  Pdo  vs  P^  for  the  linear  nodal 

is  shown  in  Figure  56  .  Inclusion  of  the  internal  resistance  term 
introduces  a  phase  shift  between  the  suction  and  discharge  pressure 
as  shown  by  the  envelope  of  operation*  The  amount  of  separation 
between  the  upswing  path  aid  the  downswing  path  was  found  to  de¬ 
pend  os  the  amplitude  of  oscillation.  If  the  same  linear  pump 
model  were  driven  by  a  suction  pressure  resulting  from  a  nonlinear 
cavitation  compliance ,  the  discharge  pressure  would  behave  as 
shewn  in  Figure  57  • 


(3)  Nonlinear  Puap  Model.  ?he  second  way  of  expressing 


puap  head  rise,  F^,  was  as  a  nonlinear  function  of  auction  pr«a- 

aure.  Data  froa  the  steady-state  head  rise  curves  were  set  into 
diode  function  generators  for  input  into  the  analog  prograa.  This 
■ethod  of  introducing  nonlinear  relationships  into  the  electrical 
analog  requires  that  such  functions  be  approxiaated  by  a  aeries 
•f  straight  line  (linear)  segments*  It  was  for  this  reason  Figures 
9%and  95  appear  as  segmented  curves.  The  nonlinear  puap  model  not 
only  produced  the  double  valued  P,  vs  P  relationship,  but  it  alao 

furnished  the  sharp  drop  in  discharge  pressure  normally  associated 
with  the  appearance  of  cavitation  within  the  punp.  (Fig.  58). 


(4)  Flight  Puap  Data.  In  ,rde?  to  evaluate  the  above 
mathematical  description  of  turbopuap  operation,  suction  and  dis¬ 
charge  pressures  recorded  during  actual  flights  of  Titan  II  vehi¬ 
cles  were  exaained. 

The  telemetered  PCM  pressure  transducer  data  were 
analysed  digitally  and  plotted  for  individual  cycles  of  oscillation. 
Beauita  of  this  comparison  of  Pd  vs  P  ,  although  cloudsd  by  noiss, 

indicate  the  puap  operation  follows  two  different  paths  under  oscil¬ 
latory  conditions  as  shown  on  Figure  59  . 


Assuming  the  pressure  transducers  have  accurately 
monitored  the  physical  situation,  a  n usher  of  observations  can  be 


ndet 

(a)  During  the  increasing  pressure  portion  of  the 
cycle,  the  relationship  between  P.  and  P  _ 

CIO  so 

is  nearly  linear*  Operation  with  decreasing 
suction  pressure  appears  to  follow  the  steady- 
state  head  rise  curve . 

(b)  The  tine  spent  between  given  points  along  the 
cycle  varies  appreciably*  The  upswing  part 
of  the  cycle  occurs  within  ~.020  sec  whereas 
the  downswing  requires  ~ .070  sec. 

(c)  The  flight  data  shows  to  be  lagging  Pfl  in¬ 
dicating  an  inherent  tine  lag  within  the  puap. 
Since  the  analog  nodal  indicated  that  P^  leads 

P  by  s  snail  anount,  it  suet  be  concluded  that 
8 

either  or  both  the  vaporisation  and  condensation 
can  not  be  considered  as  instantaneous  processes* 

V.  NONLINEAR  SYSTEM  RESPONSE 

The  electrical  analog  study  of  the  nonlinear  POGO  eodel  was 
conducted  ia  two  parts.  The  first  phase  involved  a  simplified  model 
and  served  as  a  preliminary  evaluation  of  the  proposed  nonlinear¬ 
ities  -  only  tho  oxidizer  feed  system  was  allowed  to  influence  the 
vehicle  response.  The  socond  phase  of  the  study  incorporated  the 
system  into  the  model. 

a.  Response  to  Oxidizer  System  Alone*  During  the  first 
of  the  analog  study,  it  was  founa  that  the  nonlinear  puap 
characteristics  in  conjuction  with  the  proposed  nonlinear  con¬ 
tinuity  expressions  caused  the  system  to  oscillate  far  beyond  the 
amplitudes  observed  on  the  Titan  II  flights*  For  this  reason  it 
was  necessary  to  revert  to  the  linear  puap  relationship  for  which 
limiting  was  achieved  within  the  observed  range  of  amplitudes. 

(l)  Limiting  Amplitude  vs  Resonant  Frequei,  satio. 

Figure  60  serves  to  point  out  the  values  at  which  the  gimbal 
acceleration  limited  as  a  function  of  **/*„,•  Tb*  ratio  of  th* 

natural  frequencies  is  a  relativa  indication  of  the  separation 
in  resonant  frequencies*  As  shown  in  the  figure,  the  system  can 
be  nede  to  limit  st  any  of  a  large  rang#  of  values  simply  by 
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shifting  the  resonant  frequencies  of  the  structure  or  suction  line. 

(2)  Suction  Pressure  Aaplitude  Envelope.  The  envelope 
w*  .bin  which  the  suction  pressure  oscillated  at  each  limited  con¬ 
dition  is  shown  in  Figure  6l  • 

(3)  Limiting  Amplitude  Structural  Gain.  One  of  the 
system  parameters  which  is  difficult  to  evaluate  is  the  structural 
gain.  During  this  phase  of  the  study,  the  structural  gain  was 
varied  in  an  attempt  to  determine  its  influence  on  the  total  systems 
behavior  at  the  limiting  condition.  The  results  are  shown  in 
Figure  62  • 

(4)  Pump  Truss  Flexibility.  It  was  previously  known 
that  flexibility  in  the  turhopump  supports  affected  the  stability 
of  the  system}  here  it  was  found  that  by  allowing  the  pump  to 
move  relative  to  its  structural  support  the  limiting  amplitude 
was  decreased  appreciably.  Figure  63  illustrates  the  difference 
in  response  between  the  rigid  pump  truss  and  the  flexible  pump 
truss. 


(5)  Phase  A  Conclusions,  ^he  conclusions  which  were 
drawn  at  the  end  of  the  first  phase  of  the  analog  study  were  the 
following: 


(a)  Proper  limiting  amplitudes  could  not  be  achieved 
with  the  nonlinear  pump  characteristics  because 
the  use  of  steady-state  pump  curves  to  predict 
dynamic  response  was  inaccurate.  An  attempt 
should  be  made  to  define  dynamic  pump  charac¬ 
teristics  before  making  a  final  evaluation  of 
the  models  of  the  cavitation  region. 

(b)  The  system  could  be  made  to  limit  within  the 
proper  range  of  observed  Titan  II  amplitudes 
with  each  of  the  proposed  forms  of  the  con¬ 
tinuity  expression  -  each  model  furnished  a 
different  limiting  envelope.  The  envelopes 
are  compared  in  Figure  61  and  indicate  that 
the  constant  enthalpy  process  best  achieves 
low  amplitude  limiting  va1 aes. 

b.  Response  to  Coupled  System.  The  second  phase  of  the 
analog  study  of  the  POGO  aaplitude  limiting  problem  incorporated 
the  fuel  feed  system  into  the  math  model.  The  linear  pump  charac¬ 
teristics  had  to  be  employed  to  achieve  the  desired  limiting  ampli¬ 
tudes.  Flight  data  was  evaluated  in  order  to  arrive  at  a  better 
dynamic  pump  representation. 
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oA  x  icr(gAb) 

Figure  62  Liaitiag  Qiabal  Acceleration  vs  Structural  Sain 
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(1)  Input  Parameters.  The  systems  response  was  again 
examined  at  conditions  corresponding  to  various  times  after  lift 
e?f.  Specifically  times  of  100,  110,  120,  130,  135  and  140  seconds 

jmare  investigated.  Nominal  ullage  pressures  and  propellant  con- 
'Sumption  rates  were  used.  Time  dependent  parameters  are  illustrated 
in  Figures  91*  92,  93  and  96.  Curing  the  study  some  question  was 
v  raiao*  over  which  set  of  fuel  pump  characteristics  were  appropri- 
-V-'-^  mte,  ao  pump  gains  were  parameterized  in  the  range  1.0  to  2.6. 

■a, 

(2)  General  Effects  of  Fuel  System.  The  general  patterns 
of  behavior  that  were  recorded  during  phase  A  were  again  observed 
with  the  fuel  feed  system  included  in  the  total  syatema  response. 
Because  minor  revisions  were  made  on  the  nonlinear  elements  for 

tha  second  phase  of  the  study,  it  is  not  possible  to  qualitatively 
compare  the  sole  effects  of  adding  the  fuel  ays tea. 

(3)  Pump  Gain  Sensitivity.  Because  of  lack  of  data 
describing  dynamic  pump  characteristics,  the  sensitivity  of  lim¬ 
iting  amplitude  to  this  parameter  was  investigated.  Earlier  studies 
which  deist  with  the  stability  of  the  linear  FOGO  nodal  found  that 
pump  gain  had  little  effect  at  stability.  This  was  illustrated  by 
the  two  transfer  functional 


The  expression  above  shows  a  tendenoy  for  the  pump 
gain,  H  ♦  1,  to  cancel  and  nullify  ita  effects  on  the  linear  open 
loop  gain. 

(a)  Limiting  Amplitude,  figure  64  shows  the  effect 
of  pump  gain  variation  on  the  limiting  amplitude  under  conditions 
where  the  fuel  system  response  is  strong.  It  appears  from  the 
analog  data  that  for  the  vehicle  configuration  studied,  small 
variations  in  the  open-loop  gain  produce  relatively  large  changes 
in  Halting  amplitude.  The  reason  the  limiting  amplitude  is  so 
ssnaitive  to  pump  gain  or  other  system  parameters  has  been  at¬ 
tributed  to  the  relatively  weak  nonlinearity  that  is  limiting 
the  oscillations.  For  example,  the  135  sec  configuration  that 
limited  at  H  g's  (0  to  peak)  was  oscillating  over  only  a  small 
portion  of  the  continuity  curve  a*  shown  in  Figure  65  • 
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Vigor*  65  Extent  of  Nonlinearity  -  135(sec)  United  Condition 


Were  th*  nonlinear  continuity  relationship  nore  severe ,  the  aspli- 
tude  of  oscillation  would  need  to  change  less  to  acquire  the  sane 
Mount  of  nonlinearity*  If  th*  nonlinear  continuity  expression 
were  sore  drastic,  it  nay  also  have  been  able  to  offset  th*  de¬ 
stabilising  turbopunp  nonlinearity* 

(b)  Fuel  Puap  Gain, Figure  66  shows  the  influence 
the  fuel  pssp  gain  has  on  th*  asplitude  at  which  th*  vehicle  lisits. 

c*  P Ratio*  Another  aspect  of  the  nonlinear  sodel  which 

can  be  cospared  with  flight  data  is  the  ratio  of  the  suction  pres¬ 
sure,  P  .  to  th*  structural  acceleration,  g.  Duplicating  this 
ratio  with  a  sathesatieal  sodel  isplies  th*  coupling  between  the 
structure  sad  feed  aystes  has  been  adequately  described. 

(1)  Ccsparison  of  Results .  Data  fros  some  of  the  earlier 
Titan  II  flights  is  presented  in  Figures  6?  and  68  superisposed 
on  the  flight  data  are  the  results  of  the  analog  study. 

(a)  Oxidizer  Sid*.  Although  the  P a<j/g  fros  the 
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•on  of  Analog  Data  and  Flight  Data  for  P^g 


analog  study  appeared  to  run  acaewhat  lover  than  the  flight  data, 
the  difference  between  the  analog  data  and  any  particular  flight 
in  no  More  than  variations  between  individual  flights. 

(b)  Fuel  Side.  Figure  6Q  shows  P  J%  observed  on 
the  Titan  II  flights  to  be  lower  than  that  obtaiSed  from  the  analog 
study.  Hots  that  these  ratios  are  based  on  conditions  where  Unit¬ 
ing  has  occurred.  A  closer  examination  revealed  that  ?af/g  adhered 

elosely  to  flight  data  when  the  nodel  was  oscillating  at  snail 
amplitudes  sad  was  in  the  linear  range.  As  the  aaplitudes  increased, 
and  the  nonlinearity  on  the  fuel  side  becaae  of  importance,  P  J g 
rose  tram  5*0  to  8.0  at  the  limiting  condition.  This  would  8 
indicate  that  the  fora  of  the  nonlinearity  nay  need  to  be  Modified 
on  the  fuel  aide  in  order  to  natch  thia  aspect  of  the  nonlinear 
Model  with  the  actual  vehicle. 

d.  Fuel  Feed  Syatea  Contribution.  The  effects  of  Modifying 
tbs  fuel  feed  syntes  on  the  total  systems  response  is  illustrated 
in  the  table  below.  Si  >  different  fuel  systen  representations  were 
investigated  with  reapaot  to  be  ssae  structural  and  oxidiw  systes 
configuration. 


TABLE  V  EFFECT  OF  FUEL  SYSTrH  PUMP  GAIN  AND  CONTINUITY 
EQUATION  ON  LIMITING  AMPLITUDE 


Cavitation 

Compliance 

Model 

Fuel 

Turbopuap 

Model 

Structural 
Acceleration 
Limiting  Amplitude 
(0-peak)  g*s 

Nonlinear 

Nonlinear 

2.0 

Linear 

Pusp  Gain 
-  1.0 

1.8 

Linear 

Pusp  Gain 
■  2.6 

2.7 

Nonlinear 

2.1 

Linear 

Linear 

Punp  Gain 

«  1.0 

1.9 

Linear 

Pusp  Gain 
*  2.6 

2.8 
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It  can  be  observed  from  the  table  that  use  of  a  non¬ 
linear  head  rise  aodel  on  the  fuel  side  is  de -stabilising;  use 
of  a  linear  cavitation  compliance  is  de-etabilisingj  the  amplitude 
is  not  as  sensitive  to  fuel-pump  gain  as  it  was  to  oxidiser  pump 
gain. 

5.  conclusions 

Conclusions  which  were  drawn  after  the  second  phase  of  the 
analog  study  are  the  following: 

1)  The  nonlinear  POGO  aodel  produced  a  system  that  had  a 
response  which  liaited  at  amplitudes  within  the  range 
of  values  observed  on  Titan  II  flights.  This  indicates 
that  nonlinear  phenomena  associated  with  caviatioa  in 

the  feed  systes  is  the  aost  important  system  non-linearity. 

2)  For  the  aodel  investigated,  the  closed  loop  nonlinear 

t* /atom  response  was  very  sensitive  to  the  nanner  in  which 
the  turbopuap  was  oodelad.  Examination  of  individual 
cycles  of  turbopuap  behavior  based  on  flight  data  in¬ 
dicate  possible  inadequacies  with  the  analog  dynamic 
pump  aodel.  Further  investigation  of  the  dynaaic  tur¬ 
bopuap  performance  would  be  warranted. 

3)  Because  the  cavitation  compliance  nonlinearity  was  a 
weak  nonlinearity,  small  variations  in  other  system 
parameters  produced  significant  changes  in  the  limiting 
amplitude.  Since  small  variations  in  system  parameters 
are  found  to  occur  from  flight-to-f light  one  naturally 
expects  corresponding  variation  in  POQO  limit  cycle 
amplitude  from  apparently  identical  vehicles,  "’hut  it 
would  not  be  possible  to  predict  limit  cycle  amplitudes 
on  a  given  unstable  flight. 
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SECTION  IX 


SEAL*-,  BAG  ACCUMULATOR  ANALYSIS  PROGRAM 


1.  INTRODUCTION 

The  sealed  bag  accumulator  has  been  used  with  success  on 
a  number  of  Titan  vehicle  missions.  The  accumulators  have  been 
designed  primarily  from  the  standpoint  of  the  amount  of  compli¬ 
ance  they  add  to  a  suction  system.  The  inertance, resistance, 
and  compliance  of  the  accumulator  itself  have  been  investigated 
only  over  a  limited  range  of  geometry  and  system  static  pres¬ 
sure  conditions.  As  the  amount  of  gain  or  phase  margin  pro¬ 
vided  by  tbe  present  accumulator  designs  decreases  the  dynamic 
characteristics  of  the  accumulator  will  become  more  important. 
Furthermore,  if  new  designs  become  necessary  additional  design 
data  will  be  required.  It  was  the  purpose  of  this  program 
therefore  to  investigate  the  dependence  of  accumulator  dynamic 
characteristics  on  apacific  design  parameters. 

The  first  step  of  this  effort  was  the  development  of  an 
analytical  nodel  of  the  accumulator  and  establishment  of  the 
proper  experimental  procedure  to  evaluate  the  model.  The  de¬ 
tails  of  the  analysis  phase  are  described  below. 

2.  SEALED  BAG  COMPENSATOR  ANALYSIS 
a.  Discharge  Lina  Pulsing. 

(l)  Test  configuration.  An  analysis  was  made  to 
determine  the  theoretical  response  of  the  sealed  bag  accumula¬ 
tor  test  configuration.  The  analysis  was  conducted  on  the 
eleotrioal  analog  shown  below* 
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where  i 


K  -  pump  gain 
P 

I  -  inertance,  suction  line 

I  -  inertance,  line  from  accumulator  to  pump 
ap 

I  •  inertance.  accumulator 
a 

Ig  -  inertance,  discharge  line 

C  -  compliance,  suction  line  (neglected) 

0 

C  -  compliance,  accumulator 

A 

C  *  compliance,  cavitation  of  pump,  (neglected, 
p  pump  is  not  operated  in  cavitating  region) 

-  resistance,  accumulator 

Rd  -  resistance,  discharge  line 
• 

Xp  -  pulsar  instantaneous  velocity 

-  instantaneous  suction  line  pressure 


This  configuration  is  characterized  by  the  accumulator  plaoed 
in  the  auction  line  and  excited  by  a  pulsar  located  in  the  dis¬ 
charge  line  immediately  downstream  of  the  pump. 

(2)  Results.  The  response  of  the  system  expressed 
as  suctiou  pressure  as  a  function  of  frequency  follows  the 
sketch  shown  in  Figure  70  *  The  initial  peak  or  resomanoe  is  a 
result  of  the  Interaction  of  the  suction  line  inertance,  I#, 
with  the  accumulator  compliance.  This  is  followed  by  an  anti- 
resonance  which  is  a  function  of  the  accumulator  resistance  R*r 
and  its  inertance,  I**  One  of  the  primary  objectives  of  the 
test  is  to  investigate  the  effect  of  accumulator  geometry  on 
its  resistance  and  inertance.  Obviously,  since  these  parame¬ 
ters  are  associated  with  the  anti resonance,  the  dynamic  re¬ 
sponse  will  be  too  small  to  observe.  Therefore,  the  config¬ 
uration  was  changed  so  that  the  accumulator  inertance  could 
be  interpreted  in  terms  of  effects  on  a  measurable  system 
resonance. 


U! 


frequency (epe) 


figure  70-  Typical  Response 

b.  Sealed-Bag  Compensator,  Gas  Pulsing 

(l)  Test  Configuration,  A  change  in  the  test  con¬ 
figuration  can  be  aireply  accomplished;  which  will  allow  effects 
of  accumulator  inertance  changes  to  be  evaluated  in  terms  of 
changes  to  a  system  resonant  response,  A  typical  response  is 
shown  in  the  following  figure  in  which  the  resonant  frequency, 
« jp,  is  most  sensitive  to  accumulator  inertance  and  compliance* 


m  »10.0 


frequency (cps) 

Figure  71  Typical  Qae  Puleir.g-Sealed-bag  Response 
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The  configuration  change  requires  direct  pulsing  #f 
the  gas  in  the  accimulator  bags  instead  of  the  cystesi  fluid* 

This  method  can  best  be  understood  by  referring  to  the  electri¬ 
cal  analog  diagram  below* 


Figure  72  Sealed-bag-Gae  Pulsing  Analog 


By  making  the  suction  line  as  short  as  possible, 
or  in  other  words,  by  attaching  the  accumulator  to  the  supply 
tank  outlet,  the  impedance  due  to  the  inertance,  Ia,  can  be 
made  very  small  oompared  to  the  downstream  impedance  which  in¬ 
cludes  Iap,  Id,  Bd*  With  these  impedances  in  parallel  aorose 
the  accumulator,  the  net  impedance  is  essentially  that  due  to 
Is  or  on  redrawing  the  schematic  we  hare, 


Figure  73  Simplified  Electrical  Model  For  Oas  Puleing 

The  final  result  is  a  parallel  resonant  circuit 
which  exhibits  a  high  impedance  at  resonance*  The  trans¬ 

fer  function  for  this  circuit  is  giren  by 


VA?  - 


at*.  + 12  *  \ -ja 

5(1,  *  v  *  *.  *  so- 


(61) 
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i&fcfe  A?  is  the  pulsar  voluaetrlo  displacement.  fly  keeping  I8 
pH-?*  ajpptoadLaately  «pal(  changes  in  Ia  can  readily  be  ob- 
peTVSd  la  the  aystsa  performance. 


(2)  leeults*  Epical  theoretical  values  for  systea 
yrm— ters  based  on  current  designs  give  a  response  as  shown  in 
H gor»7t*  Of  the  curves  shown,  one  represents  an  accumulator 
Laertanoe  of  *001  S#c2/in2  while  the  other  is  for  a  value  which 
I*  yOO^  greater*  These  inertance  values  cover  an  ezpeoted 
range  with  reepect  to  present  design.  Each  family  of  three 


*er**e  it  for  an  aocuaulator  resistance  which  reflects  fluid 
IWjdnf  ratios  of  0*10,  0*50,  ami  0*50.  Assuming  that  the  ac- 
aoafaater  cestriliancw,  S*,  will  he  detexalned  through  an  accurate 
jalihrationt  the  oosrpsneato r  inertanue,  I»,  can  be  aeaeured  by 
the  resoaenoe  while  the  resistance  B*,  will  be  obtained  fro*  the 
pffc  response  at  resosanoe. 
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SECTION  X 


\ SEALED  BAG  ACCIKDLATOH  TESTS 

IX  , ;  - 

'&&*■  vnxmaim 

^rS'f  ' 

9m  seeled  tag  accumulator  testa  were  performed  to  evaluate 
r  '* Ata  analytical  Moult*  of  Ssotion  IX  and  the  behavior  of  certain 
;*  Maled  tad  accumulator  dywusic  characteristics  as  a  function  of 
H  pMSfUre  and  acousulator  geometry.  Accomplishment  of  these  ob- 
,‘v,;  $W*ivee  was  attempted  through  a  test  program  of  two  phases. 

;  a  following  paragraphs  provide  a  description  of  the  test  fix- 
i  tore*  test  procedure  and  experimental  measurements  taken  during 
each  test  run.  Th*  result*  of  the  experimental  study  are  pre- 
-  tainted  in  Part  3  of  this  section. 

r*  l 

2*  TEST  PIXTOBE  AID  BSOCEDOHE 

*  . 

'  The  taste  experimental  fixture  was  the  seme  for  both  phases 

of  testing.  Figure  289  shows  the  complete  configuration  aa  an 
assembly  of  tanka,  accumulator,  suction  line,  pomp,  and  return 
or  Aibctarg*  line*  Figures  290  an4291  of  Appendix  XIV are  de¬ 
tailed  drawings  of  the  accumulator  itself  and  show  how  the  com¬ 
pliant  elements  are  looated  within  the  accumulator.  Figures 
75  through  Si  are  photographs  of  the  disasseabled  accumulator, 
including  bag*  and  retainer  sleeve.  Detailed  drawings  of  the 
rsmsi Injta  System  elements  are  given  in  Figures  290  through  292 
of  Appendix  Xiv •  Autographs  of  the  fixture  are  given  in  Fig¬ 
ures  82  sad  85  • 

a.  Sealed  Bag  Aoousulator.  The  accumulator  is  shown  as 
,  aa  element  short  coupled  to  the  supply  tank.  Aa  discussed  in 
Section  IX  above  the  short  coupling  minimises  the  inertanoe 
upstream  of  the  accumulator.  The  pulser  is  looated  dose  to 
the  aoousulator  where  it  can  effectively  exoite  the  system  by 
pulsing  th*  gas  in  the  accumulator  bags.  Figure  295  (Appendix 
XIV )  details  the  pulse  drive  assembly  which  Includes  the  vari¬ 
able  drive  motor  and  a  oam  follower  which  actuates  a  bellows- 
type  pulser  illustrated  in  Figures  292  end  293 ,  Appendix  XIV* 
figures  294  and  299  illustrate  a  valve  which  was  devised  to 
seal  off  the  inlet  side  of  the  accumulator.  Beoause  of  the 
necessity  for  short  coupling  the  accumulator  to  the  tank,  an 
ordinary  valve  could  not  be  used  at  this  location.  The  inlet 
seal  is  used  in  conjunction  with  a  ball  valve  at  the  accumula¬ 
tor  outlet  to  isolate  the  accumulator  fro a  the  remainder  of  the 
system.  This  is  done  after  static  pressures  have  been  applied 
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to  both  the  fM  in  the  ha g  and  the  tank  top  or  aye  ten  liquid* 
Ifhe  purpose  of  using  the  valves  in  this  earner  is  to  measure 
only  the  ooaplianoe  of  the  acounulator  under  operating  condi¬ 
tions  and  to  be  able  to  perform  this  operation  immediately 
preceding  a  dynamic  test. 


b*  Phase  One*  The  two  test  phases  were  integrated  in 
certain  respects.  The  first  phase  which  consisted  of  evalua¬ 
ting  the  variation  in  compliance  as  a  function  of  pressure  and 
geometry,  involved  measuring  the  change  in  volume  of  the  sealed 
accumulator  prior  to  each  of  the  second  phase  runs*  The 
method  of  measuring  volume  change  was  to  first  isolate  the 
sealed  bag  accumulator  by  meana  of  the  valves  described  above* 
The  eeoond  etep  involved  pressurising  the  bags  to  13  ASIA  and 
than  sealing  the  bags*  The  pressure  in  the  accumulator  was 
then  Increased  to  20  P5IA.  From  this  starting  point  the  pres¬ 
sure  in  the  accumulator  was  further  increased  in  5*0  PSI  in¬ 
crements  to  40*0  PPIA  and  the  volume  of  water  added  to  the 
accumulator  for  each  increment  was  measured. 


e*  Phase  Two*  The  second  phase  of  sealed  bag  accumula¬ 
tor  testing  provided  frequency  response  data  on  the  accumulator 
as  a  function  of  system  static  pressure,  per  cent  hole  area  in 
the  retainer  sleeve  and  accumulator  geometry*  The  test  pro¬ 
cedure  employed  was  very  similar  to  that  described  in  Section 
mi.  Basioally,  upon  establishing  system  steady  flow  and 
pressure  conditions  the  puleer  varldrive  system  was  actuated 
to  indues  flow  displacement  oscillations  in  the  system  with  a 
continuous  variation  in  frequency  from  2  to  JO  cycles  per  sec¬ 
ond*  The  total  time  for  the  sweep  of  frequencies  was  about  70 
seconds*  In  each  teet  the  flow  displacement  resulted  from  a 
pulsar  volumetric  displacement  of  20  cubic  centimeters*  Mea¬ 
surements  for  the  Phase  II  tests  were  pulsar  dinplacenent  and 
frequency,  flow  rate,  system  static  pressure  ana  accumulator 
oaoillatory  pressure.  Table  vi  contains  a  list  of  the  test 
runs  performed  along  with  the  pressure,  retainer  hole  area, 
and  configuration  conditions.  As  indicated  in  the  table,  four 
basio  accumulator  designs  were  examined.  In  addition,  the  bag 
retainer  hole  area  ratio  (ratio  of  hole  area  to  total  area) 
was  included  aa  a  geometric  parameter. 

The  results  of  the  sealed  bag  acounulator  test  pro- 
Iran  are  given  in  the  following  paragraphs  along  with  a  dis¬ 
cussion  of  the  conclusions  and  recommendations. 
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3.  SEALED  BAG  ACCUMULATOR  TEST  RESULTS 

/  As  desoribed  in  previous  paragraphs  the  sealed  bag  accumu¬ 
lator  tests  ears  performed  in  tvo  phases.  The  Phase  I  data 
provided  intonation  regarding  the  effects  of  accumulator  geo¬ 
metry  and  systom  static  pressure  on  accumulator  static  compli- 
aaoe,  Warn  results  of  the  Phase  II  test  along  with  the  math 
model  were  to  be  utilised  in  establishing  relationships  between 
aecaeulator  dynamic  parameters  such  as  inertance,  resistance, 
end  compliance  and  the  accumulator  design  parameters  of  geome¬ 
try,  and  system  pressure.  The  results  of  this  effort  are  con¬ 
sidered  in  the  following  paragraphs  along  with  conclusions 
derived  from  post  test  analysis. 

4.  PSASE  I  RESULTS 

la  order  to  examine  the  influence  of  accumulator  geometry 
on  compliance  four  basic  designs  were  studied.  Three  of  the 
designs  -  fuel  bag  standard  configuration,  checkout  bag  con¬ 
figuration,  and  flex  sleeve  accumulator  were  tested  with  three 
variations  in  the  retainer  sleeve  hole  area  ratio.  In  eaoh  > 
retainer  sleeves  having  hols  area  to  total  area  ratios  of  10%, 
3Q%»  and  50 %  were  examined.  The  fourth  design  was  based  on  the 
fuel  bog  configuration  but  with  the  bags  in  a  close  in  config¬ 
uration  (Figure  290  »  AppendixXIV)  and  a  10%  hole  area  ratio  in 
the  retainer*  2br  these  configurations  the  changes  in  accumu¬ 
lator  gas  volume  with  system  pressure  were  plotted}  with  the 
results  shown  in  Figures  84  through  94  •  The  bag  precharge 
pressure  in  each  case  was  15  PSIA.  In  the  figures  the  accumu¬ 
lator  compliance  at  n  given  pressure  is  related  to  the  slope  of 
the  volume  change-pressure  curve  (M)  by 

C  mi  m  (62) 

where*  »  is  the  fluid  density 

a.  Influence  of  Retainer  Hole  Area.  The  influence  of 
retainer  hole  area  on  compliance  for  a  given  design  can  be  ob¬ 
served  by  *vffi;paring  the  data  of  Figures  84^5 and  86  for  the 
fUel  bog  standard  configuration,  Fi.gures89.39  and  90  for  the 
checkout  bag  configuration,  and  Figures91,92&nd  95  for  the 
flex-sleeve  accumulator.  A  comparison  of  the  data  of  Figures 
84  through  §Q  show  that  variations  in  geometry  for  the  fuel 
hog,  both  standard  and  close  configuration,  ax.-  the  checkout 
bag  produce  only  minor  variations  in  compliance.  The  change  in 
geometry  from  these  configurations  to  that  associated  with  the 


fl«x  sleeve  accumulator  produces  a  substantial  variation  in 
compliance.  This  behavior  can  be  observed  by  comparing  the 
data  of  Figurea91t92and  93  with  that  of  Figures  84  through 
90  •  Recalling  the  relationship  between  slope  (M;  of  the  B-AT 
curve  and  compliance  (Equation  (62))  the  fl' "  sleeve  accumulator 
exhibits  more  compliance  at  a  given  system  pressure  then  do  the 
three  other  designs.  Considering  design  differences,  the  added 
compliance  of  the  flex-sleeve  configuration  can  probably  be 
associated  with  the  lower  structural  rigidity  of  the  retainer 
sleeve  coupled  with  the  bag  compliance.  With  these  conclusions 
established  the  variation  of  compliance  with  system  statio  pres¬ 
sure  for  the  fuel  bag  configurations  and  checkout  bag  can  be 
studied  independent  of  their  relative  geometries.  Ths  differ¬ 
ences  between  these  three  designs  and  the  flex-sleeve  design 
will  couple  geometry  with  pressure  effects* 

b*  Influence  of  Basic  Geometry.  Considering  Figure  84 
to  be  typical  of  the  P-  AV  curves  of  the  fuel  bag  and  checkout 
bag  configurations  it  is  apparent  that  accumulator  compliance 
decreases  with  increasing  system  pressure.  The  same  result 
applies  to  the  flex-sleeve  design  (Figures  91  through  93  )• 

This  effect  is  also  shown  in  Figure  9>4  where  the  parameter 
C/f  is  plotted  as  a  function  of  system  pressure*  Referring  to 
Figure  94  the  behavior  of  the  curves  is  not  significant  from 
the  standpoint  of  any  unique  physical  process.  The  data  osnnot 
be  related  to  either  a  gas  compression  process  or  a  structural 
deformation  process.  Considering  the  physioal  elements  con¬ 
tributing  in  the  process  a  combination  of  the  two  is  probably 
responsible  for  the  behavior  indicated, 

c.  Conclusions.  The  conclusions  drawn  from  this  phase 
of  testing  are  l)  Retainer  hole  area  ratio  does  not  signifi¬ 
cantly  affect  accumulator  static  compliance  and  2)  Accumulator 
geometry  and  system  statio  pressure  both  influence  compliance; 
however,  the  extent  of  each  is  probably  related  to  the  stiff¬ 
ness  of  the  retainer  sleeve. 

5.  PHASE  II  RESULTS 

As  was  done  in  other  phases  of  testing  the  sealed  bag  ac¬ 
cumulator  test  data  was  analyzed  in  tens  of  frequency  response* 
In  all,  four  basic  accumulator  designs  were  examined  along  with 
three  combinations  of  retainer  sleeves.  The  configuration* 
were  the  same  as  described  in  part  2  above.  Figures  97  through 
106  present  plots  of  the  experimental  P*/a?  amplitude  as  a 
function  of  frequency.  P»  ia  the  amplitude  of  the  accumulator 
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REASB  IN  ACCUMULATOR  VOLUME  -  CUBIC  INCHES 
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SEALED  BAG  ACCUMULATOR  TEST 


CHANGES  IN  ACCUMULATOR  GAS  VOLUME  WITH  PRESSURE 


Figure  93  Sealed  Bag  Accumulator  &V ,  Pressure  Relationship 


Figure  94  sealed  Bag  Accumulator  Observed  C/p 


oscillatory  pressure  and  AV  the  amplitude  of  the  puls&r  voltt* 
metric  displacement.  A  V  for  all  sealed  bag  accumulator  tests 
was  20  Cu.  CM.  Initial  examination  of  the  data  in  this  form 
indicates  that  for  the  designs  tested  both  the  frequency  cor¬ 
responding  to  maximum  response  and  the  amplitude  at  maximum 
response  increase  with  increasing  system  pressure.  This  behav¬ 
ior  suggests  that  changes  in  pressure  influence  accumulator 
compliance  primarily,  with  a  lesser  if  any  affect  on  inertance, 

a.  Correlation  with  Math  Model.  One  purpose  of  the  ac¬ 
cumulator  test  was  to  determine  if  the  observed  behavior  corre¬ 
lated  with  the  mathematical  development  described  in  Section 
IX.  The  fact  that  the  experimental  observations  do  not  match 
the  theory  is  apparent  from  a  comparison  of  the  data  of  Figure 
74  with  that  of  Figures  97  through  106.  The  lack  of  agreement 

is  found  both  in  the  frequency  range  for  maximum  response  as 
well  as  the  amplitude  of  maximum  response.  To  determine  if  the 
lack  of  agreement  was  due  to  the  values  of  parameters  used  in 
the  aralysis  of  error  in  the  basic  system  synthesis,  the  test 
data  was  replotted  aa  bode  diagrams.  Assuming  the  math  model 
to  be  correct  Figure  95  would  be  a  typical  bode  plot  of  the 
analytical  results.  Recalling  the  general  form  of  the  Pg/A  V 
transfer  function  (Equation  6:  ) 

f  -  -  (63) 

s  +  bs  +  c 

The  lollowing  information  would  be  available  from  the  test  data. 
In  Figure  V5  the  asymptote  from  the  origin  to  the  first  break 
point  (f  -  f]_)  corresponds  to  the  parameter  ^a'„  The  asymp¬ 
tote  from  the  first  break  to  the  second  break  u"  ’  fp)  yields 
the  parameter  and  finally  the  horizontal  asymptote  pro¬ 

duces  thsi  parameter  a/C  .  When  replotted,  the  test  data  did 
not  produce  the  expected  behavior.  Rather  than  approaching  a 
horizontal  asymptote  with  increr-ing  frequency  the  data  indi¬ 
cated  a  high  frequency  asymptote  with  «  slope  of  12  ED/octave 
(Reference  Figure  %  ). 

b.  Indicated  ifcdel.  This  behavior  indicates  therefore 
that  the  basic  model  synthesis  is  in  error.  If  tne  asymptotes 
indicated  by  the  test  data  are  considered  in  detail  the  ez- 
perimental  behavior  of  the  system  would  correspond  more  closely 
the  circuit  analog  given  in  Figure  io?  * 
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Figure  107  Sealed  Bag  Accumulator  Circuit  Analog 


For  this  circuit  the  E/ft  trans?*’’  function  would  be 


E  _ a _ 

q  * _ 1 _ 5 _ LZ. 

V8!  V*8! 

Equation  (64)  where  E  is  analogeous  to  Pm  and  ft 
can  be  rewritten  as 


(64) 


analogeous  to 


T 


E 

ft 


1iV  {’2 
(4  *  L2) 


n _ 2 _ 

A  *  A  « <■  1 

h  *  Y  h  *  Y 


(65) 


To  evaluate  the  experimental  parameters  for  this  function 
either  additional  test  data  would  be  required  or  the  value  of 
certain  parameters  would  have  to  be  known.  Unfortunately 
neither  set  of  information  was  available  from  the  analysis  or 
experimental  data*  The  test  data  therefore  could  be  evaluated 
on  a  qualitative  basis  only.  For  this  evaluation  the 
response  frequency  and  amplitude  were  considered  in  terms  of 
geometry  sys  am  pressure,  and  retainer  he'd  area. 
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c.  Effect  of  Retainer  r ole  Area.  The  first  relationship 
examined  was  the  effect  of  retainer  hole  area  on  the  maximum  re¬ 
sponse  frequency.  Figure  108  shows  this  effect  for  each  of  the 
accumulator  designs  with  system  pressure  as  a  parameter.  The 
data  indicates  that  for  certain  geometry  and  pressure  conditions 
a  maximal  exists  in  the  maximur  response  frequency  -  retainer 
hols  area  relationship.  Of  the  tnree  basic  designs  the  check- 

c <?t  bag  configuration  for  hole  are*  ratios  greater  than  30^ 
appears  to  have  t  higher  frequencies  associated  with  it.  The 
next  rslotinnehip  be  examined  was  that  between  the  amplitude 
of  response  and  hole  area  ratio.  Figures  109  thru  111 

show  this  effect;  respectively,  for  the  fuel  bag  standard  con¬ 
figuration,  checkout  bag  configuration  and  flex  sleeve  accumu¬ 
lator  designs.  Of  the  three  sets  of  data  only  the  fuel  bag 
standard  configuration  consistently  indicates  a  maximum  in  the 
maximum  response  amplitude.  Furthermore,  these  maximums  are 
the  highest  response  obtained  for  the  tnree  designs. 

d.  Effect  of  Pressure.  The  third  experimental  relation¬ 
ship  examined  was  the  influence  of  system  pressure  on  maximum 
response  frequency.  Figure  112  thru H4  present  plots  of  the 
experimental  data  demonstrating  this  effect.  In  each  figure 
retainer  hole  area  ratio  is  a  parameter.  A  comparison  of  the 
three  sets  of  curves  show  that  an  increase  in  system  pressure 
produces  an  increase  in  maximum  response  frequency  for  all 
three  designs  throw.^hout  the  range  of  hole  area  ratios  tested. 
The  final  qualitative  analysis  performed  on  the  data  involved 
plotting  the  accumulator  maximum  response  as  a  function  of  sys¬ 
tem  pressure  with  retainer  hole  area  ratio  as  a  parameter. 
Figures  H5 thru  117  show  these  results.  The  trends  indicated 
by  the  data  are  that  an  increasing  system  pressure  produces  an 
increase  in  amplitude  of  the  maximum  response. 

a.  Conclusions.  The  sealed  bag  accumulator  tests  pro¬ 
duced  results  which  could  be  evaluated  on  a  qualitative  basis 
in  terms  of  pressure  and  geometry  effects.  The  most  signifi¬ 
cant  test  result  was  that  the  test  data  showed  the  oath  model 
of  Section  IX  to  be  in  error  and  the  accumulator  to  be  a  more 
complicated  system  than  was  originally  anticipated.  In  order 
to  properly  evaluate  the  sealed  bag  accumulator  on  a  qualita¬ 
tive  bacis  it  wjll  be  necessary  to  revise  the  analytical  model 
along  the  lines  of  Equation  (65)  and  to  significantly  modify  the 
methods  of  testing.  Although  the  model  of  Figure  107  seems  to 
best  explain  the  experimental  data  a  number  of  parameters  in  the 
model  (I«#,  H.;»)  could  not  be  identified  in  this  program.  It  is 
net  certain  whether  the  additional  inertance  and  resistance 


terms  are  associated  with  the  accumulator;  ar  1  if  so  where  are 
they  physically  located,  or  with  the  system  external  to  the 
accumulator*  In  order  to  properly  understand  the  accumulator 
characteristics  additional  analysis  and  possibly  experimenta¬ 
tion  will  be  needed  in  the  area  of  parameter  identification. 


in 
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FREQUENCY  WITH  PER  CENT  RETAINER  HOLE  AREA 


Figure  108  Sealed  Bag  Accumulator,  Maximum  Frequency/Hole  Area 
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Figure  110  Se&led  Bag  Accumulator,  Ps/aV  with  Hole  Area 


- 1)0 

SEALED  BAG  ACCUMULATOR  BST  ~T - 

VARIATION  OF  MAXIMUM  RESPONSE  ' 

AMPT.T’nmr  inriw  WO  WW  dvmtwb  •an rv  imi 

32 

FLEX  SLEEVE  ACCUMULATOR 

— - 40  PSIA 

— 

_ 

--3U 

•—20 

_ 

rSJLA 

PSIA 

_ 

_ 

L 

_ 

_ 

_ 

_ 

_ 

■ 

1 

■ 

■ 

a 

a 

a 

a 

a 

a 

2 

a 

■ 

B 

■ 

■ 

— - 

— 

— 

a 

a 

a 

a 

a 

1 

■ 

B 

B 

« 

« 

M 

_ n 

■ 

E 

a 

a 

■ 

S 

a 

a 

a 

B 

w 

Q 

■ 

m 

a 

a 

a 

5 

> 

B 

■ 

a 

a 

a 

a 

a 

■ 

ef 

■ 

a 

a 

a 

a 

a 

a 

■ 

9 

■ 

■ 

a 

a 

B 

■ 

i 

■ 

a 

B 

B 

B 

■ 

■ 

■ 

■ 

■ 

B 

B 

a 

a 

■ 

B 

■ 

■ 

a 

a 

a 

a 

a 

■ 

a 

a 

■ 

a 

a 

a 

a 

B 

■ 

B 

a 

■ 

B 

a 

a 

a 

B 

a 

■ 

a 

B 

■ 

a 

a 

a 

a 

a 

a 

a 

a 

■ 

i 

warn 

■ 

333 

_ _ 

mm 

■■ 

S 

■ 

a 

■ 

a 

a 

a 

■ 

a 

a 

a 

■ 

a 

m 

Figure  111  Sealed  Bag  Accumulator,  P^v  with  Hole  Ar«a 


SEALED  BAG  AOCUMJLAIOR  TEST 
VARIATION  OF  MAXIMUM  RESPONSE 
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SEALED  BAG  ACCUMULATOR  TEST 
VARIATION  OF  MAXIMUM  RESPONSE 
3  Af- FREQUENCY  WITH  BAG - 


BAG  CONFIGURATION 


50%  RETAINER  HOIE  AREA 


■30%  RETAINER  HOLE  AREA 


-10%  RETAINER  HOLE  AREA 


I— 

Edom 


PRESSURE  -  FSiA 


f'igur*  11J  Sealed  Bag  Accunulatcsf ,  Maxlawa  Freque»cy/Pr««umre 
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SEALED  SAG  ACCUMULATOR  TEST 
VARIATION  OF  3&XBCK  RESPONSE 
AMPLITUDE  WHH  SAG  PREC3AI&E  PRESSURE 

FUEL  BAG  ST&MsAED  (X^FISiRATIOH 
—50%  RETAILS.  HtH£  AREA 

■ - 30%  RETAIL  SOLE  AREA 

— * - 10%  RETAINER  HQIZ  AREA 

— - 10%  RETAINER  HOLE  AREA  - 

FUEL  BAG  CLOSE  CQSFXSJRATI 


■uni 

uH 


IB— BBBl 


Figure  115  S*al»<5  Bag  Accuaulacor,  Pg/ftV  with  Pr assure 
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SECTION  XI 


CAV1TATIHG  VEHTUHI  ANALYSIS 


1.  I1TB0OTCTI0E 

Staring  the  development  of  POGO  compensation  devices  for 
the  Titan  II  vehicles  a  cavitating  venturi  was  considered  as 
a  possible  seans  of  shifting  suction  line  frequencies.  The 
sapor  phase  associated  with  the  venturi  would  produce  addition¬ 
al  compliance  in  the  system  and  thereby  lower  the  suction  sys- 
tea  natural  frequency.  Upon  examining  the  pressure  losses 
inherent  in  the  ooncept,  it  was  found  that  at  incipient  cavi¬ 
tation  the  maximum  pressure  recovery  attainable  was  on  the 
order  of  85  peroent.  To  accomplish  the  required  frequency 
shifts,  the  venturi  would  have  to  operate  in  the  range  of  well 
developed  cavitation  with  associated  pressure  losses  of  25  per¬ 
oent  or  sore.  Those  losses  were  incompatible  with  pimp  per¬ 
formance  and  tank  pressure  limitations!  hence,  the  design  was 
abondoaed.  As  a  phase  of  the  POGO  amplitude  limiting  and  control 
program  mn  alternate  to  the  pure  oavitating  venturi  was  exam¬ 
ined  for  POGO  compensation.  The  modified  concept  also  utilises 
a  oavitating  venturi;  however,  a  low  loss  tube  bypasses  the 
major  portion  of  the  miction  flow.  Figure  118  illustrates  the 
huaie  feature*  of  the  design.  Considering  that  most  of  the 
suction  flow  goes  through  the  low  Joes  bypass  tube,  the  overall 
pressure  recovery  should  be  improved  over  that  of  the  pure 
oaritatiag  venturi. 


Figure  118  Oavitating  Venturi 
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a,  Analysis.  To  determine  the  degree  of  iaprovswsnt,  an 
analysis  was  conducted  la  which  several  design*  vers  axessioed 
for  various  inlet  conditions  and  cavitation  to  bypass  flew  ra¬ 
tios.  The  basic  design  was  oriented  toward  flow  parameter* 
associated  with  the  cold-flow  test  fixture.  If  analysis  dacen- 
s t rated  the  bypass  venturi  to  be  feasible,  the  results  coaid 
be  verified  experimentally  without  additional  design  analysis* 

The  analysis  was  conducted  for  incipient  cavitation 
conditions  with  the  intent  of  extending  the  investigation  to 
fully  developed  cavitation  if,  with  the  aaxiam  idealised  pres¬ 
sure  recovery,  the  feasibility  of  the  concept  oould  be  daaon- 
strated. 


lief  erring  to  Figure  118  ,  the  flowrate  through  the 
eavitating  section  Wc  is  related  to  the  inlet  pressure  PI  and 
throat  pressure  HP  by 

-  AT. CD  ^  2gc#  (PI  -  P.r)/  1  -  Ut/AI)2  W 


where  AT  is  the  throat  area 

AX  is  the  inlet  area 

f  is  the  liquid  density 

CD  is  the  flow  coefficient  for  the  venturi 

At  incipient  cavitation  PT  ie  equal  to  the  vapor 
pressure  (PV)  of  the  fluid. 

The  recovered  pressure  (PS.) *  again  for  incipient 
cavitation,  is  given  by  (Reference  17  ) 

PS  .  Ff  +  v  (PI  -  mr)  (67) 

where  *  -  .85 

With  the  boundary  condition  of  incipient  cavitation, 
the  pressure  drop  and  thereby  the  flowrate  through  the  bypass 
section  is  established.  A  ctp-ve  of  R^yT vs  lg  was  first 
plotted  for  the  bypass  tube  (Reference  Figure  119),  R^  is  the 
Reynolds  auaber  and  f  the  friction  factor  for  the  tube.  The 
quantity  is  then  cal  .lifted  by 
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whvrs  D  is  the  tube  diameter 
Lfle  the  tube  length 
H  is  the  fluid  viscosity 

With  f  defined,  the  Reynolds  number  and  flov- 
rate  can  be  eatabliehed  from  Figure  119  •  The  streams  exiting 
the  oavitating  section  and  bypass  section  are  at  the  same 
pressure  but  have  different  velocities.  Upon  mixing  an  addi¬ 
tional  pressure  recovery  will  be  realized.  The  overall  pres¬ 
sure  recovery  is  evaluated  by  considering  the  momentum  of  the 
two  streams  at  the  tube  exit  and  at  a  station  downstream  after 
complete  mixing  has  taken  place  (Reference  18 ,  Page  221)  • 

In  addition  to  the  pressure  losses  In  the  tube  and 
Venturi  sections,  the  inlet  pressure  must  be  oorreoted  for 
approach  lueses.  These  losses  are  basically  those  associated 
with  the  contraction  of  the  fluid  stream  (Reference  18 t  Hags 

a7). 


b.  Results*  The  results  of  the  analysis  are  presented 
in  Figures  120  through  126 .  Figure  120  shows  the  variation 
of  oavitating  section  throat  area  as  a  function  of  inlet  pres¬ 
sure  and  oavitating  flow  rats  for  all  diameters  of  cypase 
tubes*  Figures  121,  122  »  *nd  125  display  the  variation  of 
bypass  flowrate  as  a  function  of  inlet  pressure  and  oaritating 
flowrate  for  bypass  tube  diameters  of  0.5  in.,  O.625  in.  and 
0*75  in.,  respectively.  The  pressure  recovery  rasio  ie  given 
in  Figures  124 »  125  *  and  126  • 

e*  Conclusions*  The  feasibility  of  the  bypass  venturi 
as  a  FOGO  compensation  devios  is  established  by  Figures  124 
12 5f  and  126  •  Since  the  maximum  attainable  recovery  pressure 
ratio  is  no  more  than  90  peroent,  which  is  at  the  limit  of 
acceptable  losses ,  fully  developed  cavitation  will  easily 
increase  the  pressure  losses  to  unacceptable  levels.  Although 
the  analysis  just  described  applies  to  a  subscale  system  com¬ 
pared  to  the  flight  system,  the  Reynolds  number  and  cavitation 
to  bypass  flow  ratios  would  be  oompar&ble.  The  systems  there¬ 
fore  will  also  be  similar  in  terms  or  recovery  pressure  ratios. 
Henoe,  from  the  analysis  results,  it  was  concluded  that  the 
bypass  oavitating  venturi  would  not  be  a  practical  POGO  com¬ 
pensation  device.  Considering  the  range  of  flight  inlet  pres¬ 
sures  over  which  such  a  fix  would  have  to  operate,  the  per¬ 
formance  would  be  compromised  even  further. 
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